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ABSTRACT 

Studies  of  new  array  data  are  reported,  with  particular  emphasis 
on  the  array  NORSAR.  Models  for  approach  angle  anomalies  at 
LASA  are  presented.  Further  evidence  for  deep  mantle  lateral 
heterogeneity  is  given.  Detailed  analyses  have  been  made  of 
earthquakes  in  Eurasia.  A  microbarographic  noise  study  has 
been  completed.  Continued  progress  is  reported  in  the  bulletin 
production  from  the  month  of  intensive  seismic  analysis. 
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SUMMARY 


This  is  the  eighteenth  Semiannual  Technical  Summary  report  of  the  Lincoln  Laboratory 
Seismic  Discrimination  Group.  Our  work  is  aimed  at  the  solution  of  seismological  problems 
concerned  with  the  detection,  location  and  identification  of  seismic  events  in  order  to  provide 
technological  background  that  would  be  needed  in  discussions  of  a  test-ban  treaty.  We  have 
rapid  access  to  digital  data  from  the  large  arrays  LASA,  NORSAR  and  ALPA  in  addition  to  film 
data  from  the  World-Wide  Network. 

NORSAR  has  now  been  in  full  operation  for  more  than  a  year  and  we  have  devoted  consider¬ 
able  effort  to  studying  its  short  period  P-wave  characteristics.  These  include  a  greater  degree 
of  incoherence  and  amplitude  fluctuation  from  site  to  site  than  is  observed  at  LASA.  Our  first 
concern  has  been  to  accumulate  a  sufficient  body  of  data  to  see  whether  there  are  any  clear 
patterns  to  these  fluctuations.  It  seems  at  first  sight  that  there  are  not.  Evidence  is  presented 
that  sources  relatively  close  to  one  another  have  substantially  different  amplitude  patterns  across 
NORSAR.  We  foresee  that  an  important  effort  in  the  near  future  will  be  an  attempt  to  determine 
the  extent  to  which  fluctuations  can  be  modeled  in  a  geophysically  satisfying  way.  Many  seismic 
source  regions  are  within  30°  of  NORSAR,  and  accurate  location  of  these  requires  a  good  under¬ 
standing  of  the  dT/dA  curve  in  this  range.  This  is  the  distance  range,  however,  for  which  we 
expect  most  variability  owing  to  crustal  and  upper  mantle  lateral  heterogeneity.  We  report  newly 
measured  dT/dA  values  in  this  range  and  remark  on  the  similarity  of  the  curve  to  that  for  Western 
North  America. 

Work  is  in  progress  on  the  time  delays  within  subarrays  at  NORSAR.  This  is  to  consider  the 
problem  of  whether  beaming  within  individual  subarrays  needs  to  be  more  sophisticated  than  a 
simple  plane  wave  approach. 

Collaboration  with  the  scientific  staff  of  NORSAR  continues  at  a  high  level.  The  most  recent 
visitor  to  NORSAR,  Mr.  Sheppard,  is  now  in  the  process  of  compiling  an  extensive  set  of  time 
corrections  for  subarrays  of  NORSAR.  This  will  be  important  for  two  reasons:  (1)  it  will  con¬ 
tribute  to  the  refinement  of  locations  with  NORSAR  and  (2)  it  will  enable  us  to  search  for  the 
sort  of  propagation  anomalies  which  were  found  at  LASA  and  associated  with  heterogeneity  in  the 
middle  and  deep  mantle. 

If  anisotropy  in  elastic  properties  exists  in  the  crust  and  upper  mantle  -  a  frequent  conjec¬ 
ture  as  yet  only  sparsely  documented  -  it  might  be  possible  to  use  an  array  to  look  at  surface 
wave  velocity  variability  as  a  function  of  azimuth.  A  search  has  been  made  for  anisotropy  both 
at  IJ\.SA  and  NORSAR,  but  the  results  so  far  emphasize  the  very  serious  problem  of  separating 
anisotropy  from  heterogeneity  both  within  the  array  and  beyond  it.  Frequency  dependence  of  the 
approach  angle  of  Rayleigh  waves  to  LASA  is  associated  with  continental  complexity  to  the  west 
of  the  array. 

An  attempt  is  reported  to  invert’the  large  set  of  (dT/dA,  Azimuth)  data  collected  from  indi¬ 
vidual  subarrays  at  LASA  by  placing  a  highly  corrugated  interface  at  depth.  This  surface  cannot 
be  ascribed  to  any  particular  depth  beneath  the  array,  but  if  it  were  associated  with  the  Moho  it 
could  be  mapped  out  in  an  approximate  form. 

An  alternative  approach  to  the  problem  of  interpreting  this  data  set  is  to  treat  the  scatter  in 
the  data  as  being  indicative  of  a  medium  with  scattering  properties  describable  by  a  Chernov- 
type  model.  This  technique  is  demonstrated  and  it  is  shown  how  certain  parameters  of  a  random 
medium  could  be  inferred  from  the  data.  Further,  the  impact  of  a  random  element  on  seismic 
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array  techniques  for  event  location  is  discussed.  Some  events  at  LASA  appear  more  coherent 
than  others,  and  in  one  particular  case  of  an  incoherent  event  we  have  been  able  to  identify  two 
separate  P- waves  arriving  simultaneously  at  the  array  from  substantially  divergent  directions. 

We  believe  that  this  multipathing  evidence,  rather  similar  to  that  obtained  some  time  ago  for 
Rayleigh  waves,  may  lead  us  to  further  interpretation  of  the  middle  mantle  as  a  laterally  heter¬ 
ogeneous  medium.  Work  has  advanced  on  an  extensive  travel  time  survey  of  the  mantle  using 
deep  focus  earthquakes.  It  is  possible  to  map  some  regions  of  the  mantle,  particularly 
near  the  core  boundary,  which  have  seismic  velocities  differing  significantly  from  the  mean. 

A  study  has  been  completed  of  the  contribution  of  S-waves  to  the  coda  of  a  P  arrival.  It  has 
been  possible  to  infer  a  structure  (probably  at  the  sediment-hard  rock  interface)  which  could 
generate  such  P-S  converted  energy. 

A  survey  of  focal  mechanisms  of  Eurasian  earthquakes  is  under  way.  This,  it  is  hoped, 
will  give  us  a  clearer  understanding  of  the  events  against  which  we  have  to  discriminate  in  that 
region.  It  may  also  assist  us  in  the  explanation  of  so-called  "anomalous"  events  which  fail  the 
M^m^  discriminant.  Our  initial  results  highlight  the  extraordinary  diversity  of  mechanisms 
away  from  well-defined  plate  boundaries.  A  special  study  has  been  completed  of  a  large 
Persian  earthquake  for  which  parameters  such  as  stress  drop,  dimensions  and  slip  could  be 
inferred.  Also  reported  are  the  statistics  of  a  major  earthquake  swarm  for  which  many  singular 
properties  have  been  noted.  Comparison  is  made  of  the  energy  release  time  profile  with  recent 
work  of  Knopoff  on  stochastic  seismicity.  A  deep  focus  earthquake  in  the  Hindu  Kush  region  has 
been  examined  for  its  rather  unusual  suppression  of  the  depth  phase  pP.  Without  this  phase, 
there  would  be  grounds  for  suspicion  over  this  event,  but  the  long  period  phase  sP  helps  to  con¬ 
firm  the  network-determined  depth.  Further  examination  of  a  Siberian  earthquake  by  means  of 
time-domain  analysis  of  WWSSN  seismograms  is  described. 

Work  on  the  correlation  between  microbarographic  and  long  period  seismic  noise  has  been 
completed.  It  is  possible  to  suppress  seismic  noise  by  several  dBs  with  cross-correlation 
techniques,  but  doing  this  on  a  routine  basis  requires  continued  re-evaluation  of  cross  correla¬ 
tion  functions  as  the  coupling  between  atmosphere  and  ground  is  dependent  on  wind  speed  and 
direction.  Burial  of  instruments  at  200  m  may  be  in  the  long  run  a  more  economical  way  of 
suppressing  this  noise.  The  question  of  optimal  monitoring  of  surface  waves  at  distances  of 
less  than  20  degrees  is  under  consideration.  We  report  signal  and  noise  characteristics  at 
LASA  when  explosions  at  NTS  are  observed.  At  this  distance  (12°)  the  signal-to-noise  ratio 
appears  relatively  constant  from  12  to  20  seconds. 

Work  on  the  International  Seismic  Month  continues  to  progress.  At  the  moment,  we  are 
examining  in  great  detail  the  detection  logs  of  the  LASA  and  NORSAR  arrays  in  conjunction  with 
all  other  single  station  data  in  order  to  look  for  credible  events  to  add  to  our  present  bulletin. 
Substantial  interactive  computer  aid  is  now  provided  by  a  new  system  called  DADS  modeled  on 
our  experience  with  the  display  console. 

Finally,  we  report  two  more  general  seismological  studies  performed  in  the  laboratory  in 
the  last  6  months.  A  study  of  core  reflections  PcP  has  shown  the  constraints  within  which  core/ 
mantle  boundary  models  must  lie.  Work  on  free  oscillation  has  shown  the  value  of  using  the 
focal  mechanism  of  an  earthquake  in  extracting  spectral  peaks.  Comparison  of  peaks  from  the 
Alaskan  and  Colombian  quakes  shows  systematic  and  significant  divergences  which  probably 
have  to  be  associated  with  mantle  heterogeneity. 


D.  Davies 
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ALPA 

Alaskan  Long- Period  Array 

CGS 

United  States  Coast  and  Geodetic  Survey 
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Data  Analysis  and  Display  System 
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Large  Aperture  Seismic  Array,  Montana 
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International  Seismic  Month 

NORSAR  Norwegian  Seismic  Array 
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National  Ocean  Survey 
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Nevada  Test  Site 

SATS 

Semiannual  Technical  Summary 

WWSSN 

World-Wide  Standard  Seismographic  Network 
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SEISMIC  DISCRIMINATION 


I.  ARRAY  STUDIES 

A.  SHORT- PERIOD  AMPLITUDE  AND  WAVEFORM  STUDIES  AT  NORSAR 

Variations  of  short  period  amplitudes  and  waveforms  across  NORSAR  are  being  studied  as 
part  of  the  NORSAR  evaluation  program.  Although  the  data  are  complex  and  often  inconsistent, 
some  general  conclusions  can  be  made. 

Examination  of  individual  sensors  within  subarrays  shows  that  subarray  amplitudes  are  con¬ 
trolled  by  the  amount  of  complex  scattering  of  incoming  waves  by  the  crust  under  NORSAR.  Sub¬ 
arrays  which  have  high  amplitudes  have  coherent  waveforms  within  each  subarray.  Low  ampli¬ 
tude  subarrays  on  the  other  hand  often  have  incoherent  waveforms  with  complex  multiple  arrivals 
in  the  waveform  codas.  Interestingly  enough,  on  both  low  and  high  amplitude  subarrays,  the 
center  seismometer  trace  is  often  very  coherent  with  the  steered  subarray  beam.  This  suggests 
that  center  sensors  can  be  used  in  place  of  subarray  beams  for  routine  seismic  analysis  of  events, 
provided  the  signal-to-noise  ratio  (SNR)  is  adequate. 


TABLE  1-1 

EVENTS  USED  IN  NORSAR  SHORT- PERIOD  STUDY 


Event 

Date 

Origin 

Time 

(GMT) 

Region 

Azimuth 

(deg) 

Distance 

(deg) 

Depth 

(km) 

Magnitude 

(mb) 

1 

2-21-72 

19:34:50 

Alaska  Peninsula 

-6.2 

63.3 

60D 

5.7 

2 

8-16-71 

04:58:00 

Szechwan  Province 

71.8 

67.1 

33N 

5.5 

3 

6-30-71 

03:56:57 

Eastern  Kazakh 

73.2 

39.0 

0G 

5.4 

4 

5-25-71 

04:02:58 

Eastern  Kazakh 

74.0 

38.7 

0G 

5.2 

5 

6-15-71 

22:15:50 

Kirgiz-  Sinkiang 
border  region 

82.0 

45.3 

33N 

5.1 

6 

8-4-71 

01:59:03 

Hindu  Kush 

93.7 

44.8 

206 

5.0 

7 

6-26-71 

22:23:29 

Afghanistan-  USSR 
border  region 

94.7 

44.9 

127 

5.0 

Here  we  give  an  example  of  these  results  for  a  high  frequency,  impulsive  event,  a  presumed 
explosion  in  Eastern  Kazakh,  on  25  May  1971,  with  origin  time  04:02:58  GMT.  Other  data  for 
this  event  are  given  in  Table  1-1.  P-wave  beams  at  each  subarray  were  formed  by  aligning  the 
individual  sensors.  The  peak-to-trough  beam  amplitudes  (A/T)  were  measured  and  normalized 
by  their  geometric  mean  across  NORSAR.  These  amplitudes  are  plotted  in  dB  relative  to  their 
geometric  mean  in  Fig.  1-1.  Center  sensor  amplitudes  are  also  plotted  relative  to  their  geo¬ 
metric  mean.  In  spite  of  a  variation  of  20  dB,  that  is,  an  amplitude  fluctuation  of  10,  across  the 
whole  NORSAR,  the  center  sensors  are  nearly  always  within  1  or  2  dB  of  the  beam  amplitudes. 
Also  shown  in  Fig.  1-1  is  the  SNR  at  2  Hz  measured  for  each  subarray  beam.  The  fact  that  the 
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SNR  pattern  closely  follows  the  signal  amplitude  pattern  demonstrates  that  true  P-wave  fluctua¬ 
tions  are  being  measured,  not  problems  with  seismometer  gain  or  coupling  with  the  ground.  For 
this  event,  subarrays  which  have  high  amplitudes,  such  as  4C  and  6C,  have  very  coherent  signals 
within  each  subarray.  Figure  1-2  shows  the  individual  sensors  and  subarray  beams  for  4C  and 
6C,  as  examples.  Figure  1-3  shows  similar  traces  of  two,  low-amplitude  subarrays,  7B  and  8C. 
Within  each  subarray  the  relative  gains  of  single  sensors  and  beams  are  preserved. 

There  are  significant  differences  between  the  high  and  low  amplitude  subarrays  shown  here. 
The  individual  traces  of  4C  and  6C  in  Fig.  1-2  are  very  coherent  with  rapidly  decaying  codas  in 
the  5.0- sec  window  indicated.  The  first  arrivals  are  impulsive  and  the  largest  in  the  window. 

The  first  arrivals  on  subarrays  7B  and  8C  in  Fig.  1-3  are  weak  and  the  codas  very  complex  show¬ 
ing  distinct  multiple  arrivals  with  different  phase  velocities  from  the  first  motions.  This  suggests 
that  the  waves  have  been  backscattered  by  complex  inhomogeneities  under  the  array,  the  energy 
being  distributed  among  complex  arrivals  in  the  coda.  The  center  sensors,  7B0  and  6C0,  corre¬ 
late  quite  well  with  the  corresponding  subarray  beams.  Correlation  coefficients  were  computed 
between  individual  sensors  and  the  subarray  beams  using  5  sec  of  data.  Over  the  frequency  band 
1.5  to  3.5  Hz,  7B0  and  the  7B  beam  have  a  correlation  coefficient  of  0.89,  whereas  8C0  and  the 
8C  beam  have  a  correlation  coefficient  of  0.87.  Other  sensors  have  correlation  coefficients  as 
low  as  0.31  with  the  subarray  beam.  This  indicates  that  the  arrivals  in  the  coda  are  correlated 
over  at  least  the  spacing  between  subarray  sensors  (~3  km),  and  that,  qualitatively,  center  sen¬ 
sors  can  be  taken  as  approximations  to  the  subarray  beam. 

Under  this  assumption,  amplitude  fluctuations  for  seven  events  were  measured  using  the 
center  sensors.  These  events  are  listed  clockwise  by  azimuth  in  Table  1-1.  Except  for  event  1, 
the  azimuths  sample  a  narrow  sector  which  includes  the  Eastern  Kazakh  test  site  in  the  USSR. 
Amplitudes  were  measured  and  divided  by  the  geometric  mean  of  each  event  over  the  whole  array. 
Steered  beams  across  NORSAR  are  displayed  in  Fig.  1-4,  to  indicate  the  variety  of  complexity  and 
period  of  the  events. 

Figures  1-5  through  -11  are  contour  maps  of  amplitude  fluctuation  across  NORSAR  for  the 
seven  events.  The  amplitudes  are  contoured  in  dB  above  the  geometric  mean  of  each  event,  with 
a  contour  interval  of  2  dB.  Superposed  on  each  map  is  an  arrow  which  extends  from  the  1A  sub¬ 
array  along  the  event  azimuth. 

The  most  interesting  feature  of  these  maps  is  the  pocket  of  low  amplitudes  which  appears 
on  all  but  two  of  the  contour  maps.  This  pocket  moves  about  over  NORSAR  for  different  events, 
which  suggests  that  it  may  be  caused  by  inhomogeneities  deep  in  the  crust  or  upper  mantle,  at 
least  several  times  deeper  than  the  distance  between  subarrays.  However,  if  we  put  a  hypotheti¬ 
cal  obstacle  at  depth  to  block  out  and  scatter  the  incoming  waves,  the  movement  of  the  low  am¬ 
plitude  pattern  with  azimuth  does  not  seem  to  be  consistent  for  these  events.  Another  difficulty 
is  the  consistency  of  amplitude  patterns  for  events  with  equal  azimuths.  Figures  I- 10  and  -11 
show  contour  maps  for  two  Hindu  Kush  events  which  differ  in  depth  by  more  than  100  km.  The 
amplitude  patterns  are  quite  similar,  event  7  in  Fig.  I- 11  showing  a  slightly  deeper  pocket  than 
event  6  in  Fig.  I- 10.  Each  event  shows  a  consistent  increase  in  amplitude  of  at  least  15  dB  from 
subarrays  1A  to  3C,  in  the  northeasterly  direction.  Since  the  two  events  have  different  depths, 
the  amplitude  patterns  must  be  associated  with  structure  under  NORSAR,  rather  than  source 
variations. 

1 

Lacoss  analyzed  a  large  population  of  beamed  subarray  amplitudes  generated  by  the  Event 
Processor  at  NORSAR.  One  important  conclusion  of  the  study  was  that  the  large  amplitude  fluc¬ 
tuations  across  NORSAR  are  often,  but  not  always,  repeatable  for  events  from  the  same  region. 
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As  a  counter  example,  the  two  presumed  explosions,  events  3  and  4,  in  Figs.  1-7  and  -8, 
respectively,  do  not  give  similar  amplitude  fluctuations,  in  spite  of  being  very  close  in  azimuth 
and  distance  from  NORSAR.  As  shown  in  Fig.  1-4,  the  beam  for  event  4  is  high  frequency  and 
impulsive,  and  similar  in  frequency  content  to  the  Hindu  Kush  events.  Except  for  subarrays 
4C  and  6C,  which  have  high  amplitudes,  the  contour  map  for  event  4  is  quite  similar  to  those  for 
events  6  and  7.  Event  3,  however,  seems  to  have  a  complex,  low  frequency  component  in  the 
source  function,  as  shown  in  Fig.  1-4.  The  amplitudes  were  measured  on  the  second  pulse, 
which  seems  to  contain  interfering  pulses  on  single  sensors.  Clearly,  amplitude  patterns  from 
such  an  event  cannot  be  attributed  only  to  structure  under  NORSAR. 

Event  2  is  a  Chinese  earthquake  which  has  an  azimuth  close  to  that  of  events  3  and  4,  but  is 
30°  further  away  from  NORSAR.  The  contour  map  for  this  event  is  shown  in  Fig.  1-6.  In  this 
case,  the  low  amplitude  region  is  displaced  to  the  north  relative  to  the  pattern  of  event  4  in 
Fig.  1-8.  This  is  not  easily  explained  by  the  small  difference  in  azimuth  between  the  two  events. 

These  examples  illustrate  the  complexity  of  the  amplitude  fluctuations  at  NORSAR.  Con¬ 
siderably  more  data  must  be  examined  before  the  causes  of  the  fluctuations  can  be  systematically 
described  or  explained.  c  w  Frasier 

B.  ^  MEASUREMENTS  USING  NORSAR  IN  THE  DISTANCE  RANGE  FROM  0  to  30° 


In  order  to  use  NORSAR  to  detect  and  locate  events  in  the  western  Soviet  Union,  it  is  neces¬ 
sary  to  use  P-waves  recorded  at  distances  less  than  30°,  which  are  subject  to  various  complica¬ 
tions  caused  by  the  fact  that  the  travel  time  curve  is  multi-valued  and  regionally  varying.  Fail¬ 
ure  to  account  for  these  effects  will  lead  to  large  location  errors  and  could  even  cause  one  event 
to  be  mistaken  for  several  events  in  different  places.  Therefore,  we  have  begun  "calibrating" 
the  upper  mantle  in  the  vicinity  of  NORSAR  by  making  dT/dA  measurements  for  P  waves  at  dis¬ 
tances  less  than  30°.  Figure  1-12  shows  the  data  obtained  so  far,  together  with  the  data  meas- 
2 

ured  by  Johnson  using  the  extended  Tonto  Forest  array  in  Arizona,  and  a  theoretical  curve 
appropriate  to  western  North  America.  No  attempt  has  been  made  to  apply  any  corrections  for 
the  effects  of  the  structure  beneath  the  array.  This  will  be  done,  and  azimuth  as  well  as  dT/dA 
anomalies  will  be  considered,  after  more  measurements  have  been  made. 

Already,  it  is  apparent  that  the  structure  of  the  upper  mantle  is  not  grossly  different  in 
Europe  than  it  is  in  western  North  America.  The  only  obvious  anomalies  in  the  NORSAR' data 
are  low  dT/dA  values  between  1000  and  2000  km,  probably  indicative  of  higher  velocities  in  the 
upper  250  km  in  Scandinavia,  and  a  few  scattered  points  between  9  and  10  sec/deg  (where  the 
North  American  data  are  sparse).  When  a  greater  number  of  data  have  been  analyzed  at  greater 
length,  smaller  scale  structural  details  may  emerge. 

B.  R.  Julian 


C.  TIME  DELAYS  WITHIN  NORSAR  SUBARRAYS 

There  are  at  least  two  important  respects  in  which  short  period  P-wave  data  recorded  at 
NORSAR  differ  from  those  recorded  at  LASA:  much  higher  frequencies  are  often  observed  and 
events  of  interest  are  often  much  closer  to  the  array.  NORSAR's  "window"  for  high  frequencies 
is  potentially  one  of  its  most  valuable  characteristics,  but  in  order  to  take  advantage  of  it,  time 
corrections  for  the  individual  sensors  must  be  accurately  known.  Otherwise,  the  high  frequen¬ 
cies  may  be  degraded  by  the  process  of  beamforming.  This  is  obviously  true  for  the  array  as 
a  whole,  but  it  is  less  clear  whether  delays  within  individual  subarrays  also  need  to  be  considered. 
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Figure  1-13  illustrates  that  they  do,  at  least  in  some  cases.  It  shows  beams  of  the  P  wave 
from  a  deep  Hindu  Kush  earthquake  formed  using  a  single  subarray  (2B),  using  both  hand-picked 
time  delays  and  delays  calculated  assuming  a  plane  wavefront  with  the  velocity  and  azimuth 
determined  by  fitting  the  arrival  times  at  all  the  subarrays.  The  second  beam  has  clearly  dis¬ 
torted  the  waveform  to  an  extent  which  could  hamper  its  usefulness  in  source  mechanism  studies. 

In  order  to  determine  the  seriousness  of  this  problem,  time  delays  have  been  measured  at 
all  NORSAR  sensors  for  12  events  with  clearly  recorded,  impulsive  P  waves.  The  distribution 
of  these  signals  in  slowness  space  (Fig.  1-14)  was  chosen  to  be  as  wide  as  possible,  so  that  the 
variations  in  the  delays  would  be  extreme. 

Figure  1-15  shows  the  expected  spectral  loss  caused  by  beamforming  as  a  function  of  the 
root-mean- square  (rms)  time  error.  It  illustrates,  for  example,  that  if  we  wish  the  degradation 
of  the  signal  to  be  less  than  5  dB  at  all  frequencies  less  than  5  Hz,  then  the  rms  time  error  must 
be  0.03  sec,  or  less.  (This  calculation  is  based  upon  the  assumption  that  the  signals  from  the 
individual  sensors  are  identical.  In  reality,  incoherence  would  unavoidably  degrade  real  signals 
even  further.)  Since  short- period  NORSAR  data  are  digitized  at  0.05  sec  intervals,  under  the 
best  imaginable  circumstances  timing  errors  would  be  approximately  uniformly  distributed  be¬ 
tween -0.025  and  +0.025  sec  and  the  rms  error  would  be  slightly  less  than  0.015  sec. 

The  distribution  of  the  rms  errors  determined  from  the  measured  times  is  given  by  the 
histogram  in  Fig.  I-l6(a).  Roughly  half  the  sensors  have  rms  errors  greater  than  0.04  sec,  so 
it  is  clear  that  some  sort  of  time  corrections  are  needed  if  full  advantage  is  to  be  taken  of  high 
frequency  signals.  Figure  I- 16(b)  shows  the  improvement  in  distribution  which  can  be  obtained 
by  using  a  constant  time  correction  (independent  of  velocity  and  azimuth)  for  each  sensor.  Still 
nearly  40  percent  of  the  sensors  have  rms  errors  exceeding  0.04  sec.  The  improvement  which 
may  be  gained  with  distance  and  azimuth  dependent  corrections  is  presently  being  investigated. 

B.  R.  Julian 


D.  EVALUATION  OF  NORSAR  TIME  ANOMALIES 

A  data  base  for  the  study  of  the  NORSAR  P-wave  time  anomalies  has  been  collected.  These 
data  were  collected  from  plotted  outputs  of  the  event  processor  runs  at  NORSAR  between  mid 
1971  and  mid  1972.  The  initial  selected  population  consisted  of  all  events  reported,  both  on  the 
PDE  cards  and  on  the  NORSAR  station  bulletins,  with  a  quality  rating  of  1  or  2.  From  this 
group  of  over  1600  events,  a  final  data  base  of  721  events  was  chosen  as  meeting  recording 
quality  for  reliable  delays.  Each  event  in  the  initial  group  was  examined  for  proper  delays,  as 
recorded  in  the  subarray  beam  alignment  on  the  output  plots  of  the  event  processor.  If  any  er¬ 
rors  in  alignment  were  observable  on  the  plots,  then  these  data  were  corrected  before  punching 
the  data  on  cards.  The  data  recorded  on  cards  include  all  event  information  and  array  delays. 

In  this  manner,  a  collection  of  events  recorded  at  NORSAR  over  a  13  month  period  was  examined, 
corrected  and  punched  on  cards  that  now  consist  of  the  data  base  for  NORSAR  subarray  time 
anomaly  study. 

Since  60  percent  of  these  events  in  the  final  population  need  some  corrections  made  to  the 
delays  determined  by  the  event  processor,  the  entire  data  base  has  been  entered  into  the  DADS 
system  reported  in  Sec.  V-B.  This  system  has  features  which  will  facilitate  the  updating  of 
corrections  required  for  the  population.  Additional  events  will  be  added  as  soon  as  sufficient 
numbers  of  NORSAR  event  recordings  are  received. 

R.  M.  Sheppard 
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E.  CHARACTERISTICS  OF  10  TO  40  SECOND  RAYLEIGH  WAVES 

AT  LASA  AND  NORSAR 

The  frequency- wavenumber  spectra  of  vertical,  long-period  LASA  and  NORSAR  recordings 

of  Rayleigh  waves  have  been  computed  to  determine  the  propagation  effects  due  to  crustal  struc- 
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ture  along  the  paths  to  and  beneath  the  arrays.  High  resolution  analysis  was  used  to  obtain 
the  greater  angular  resolution  and  accuracy  that  the  method  provides  over  standard  techniques. 
Phase  velocity  curves  for  both  arrays  were  determined.  Lateral  refraction  as  a  function  of 
wavelength  for  the  path  NTS  to  LAS  A  was  examined.  Velocity  anisotropy  of  the  crust  beneath 
NORSAR  was  investigated.  Multipaths  from  various  epicenters  to  NORSAR  were  also  determined. 

The  data  base  for  LASA  is  presented  in  Table  1-2.  It  is  composed  of  seven  NTS  explosions 
and  two  regional  earthquakes.  More  detailed  information  on  the  explosions  may  be  found  in 
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Springer  and  Kinnaman.  The  explosions  were  chosen  because  they  have  virtually  identical  epi¬ 
centers,  thus  propagation  effects  between  NTS  and  LASA  could  be  detailed.  The  two  earthquakes 
were  added  to  complete  the  measurement  of  the  local  phase  velocity  curve.  The  data  base  for 
NORSAR  is  presented  in  Table  1-3.  These  events  were  picked  because  they  are  well  distributed 
in  azimuth,  thus  the  effect  of  the  structure  beneath  this  array  could  be  determined. 

For  each  event,  a  high  resolution  horizontal  wavenumber  power  spectrum  at  a  set  of  periods 
spanning  the  10-  to  40-sec  band  was  computed  and  contoured.  The  wavenumber  where  the  peak 
(or  peaks  in  cases  of  multipathing)  of  the  spectrum  occurred  was  determined,  yielding  the  azimuth 
and  phase  velocity  of  the  particular  period  component  of  the  Rayleigh  wave.  Care  was  taken  in 
excluding  windowed  energy  from  neighboring  periods  with  higher  power  making  use  of  the  sensi¬ 
tivity  function  of  the  filter  used  to  isolate  the  period  under  study.  The  peak  wavenumber  was 
chosen  to  be  at  the  geometrical  center  of  the  2-dB  contour.  The  local  maximum  point  in  the 
wavenumber  spectrum  always  fell  within  one  differential  wavenumber  unit  of  the  chosen  point. 

Table  1-4  presents  the  phase  velocity  data  for  LASA  and  NORSAR.  Each  point  represents 
the  average  of  9  values  for  LASA  and  6  values  for  NORSAR.  Accuracy  is  a  function  of  period, 
being  about  ±0.05  for  short  periods  and  monotonically  increasing  to  ±0.1  for  longer  periods. 

Figure  1-17  is  a  plot  of  the  LASA  values  superimposed  on  average  curves  for  mid-continent  and 
basin- range  structures.5  The  jump  in  the  curve  near  20  sec  indicates  a  rise  in  the  shear  veloc¬ 
ity  in  the  middle  crust.  Otherwise,  the  curve  tends  to  average  the  basin-range  and  mid-continent 
curves.  Figure  1-18  gives  the  NORSAR  curve  plotted  with  an  average  shield  curve5  and  a  curve 
for  Finland.^  The  most  distinctive  feature  is  the  dip  in  phase  velocity  near  35  sec.  Waves  with 
this  period  sense  the  uppermost  mantle  so  that  lower  than  normal  shield  shear  velocities  may 
exist  there.  The  NORSAR  area  does  deviate  from  most  shield  areas  in  that  a  Permian  pluton 
fills  the  large  Oslo  graben  on  which  the  array  is  partially  located.  In  general,  these  data  deter¬ 
mined  for  NORSAR  and  LASA,  represent  the  highly  accurate  measurements  that  can  be  made 
over  relatively  short  paths  (LASA  spans  approximately  one  wavelength  of  a  40- sec  Rayleigh  wave) 
by  combining  direct  digital  recording  of  an  array  of  identical  instruments  with  numerical 
techniques. 

Julian  has  demonstrated  via  ray-theoretical  calculations  how  surface  waves  will  be  re¬ 
fracted  by  lateral  inhomogeneities  in  the  crust.  The  path  between  NTS  and  LASA  is  geologically 
complicated  in  that  it  passes  near  to  or  within  the  Great  Nevada  Basin,  the  Snake  River  Down- 
warp  and  the  Colorado  Plateau.  Each  of  these  geologic  provinces  has  different  velocity  struc¬ 
tures7  so  that  the  transition  zones  may  act  as  refractors  to  Rayleigh  waves.  Further,  any  de¬ 
flection  of  a  Rayleigh-wave  wavefront  will  be  frequency  dependent,  since  the  velocity  differences 
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TABLE  1-2 

DATA  BASE  FOR  LASA  RAYLEIGH- WAVE  STUDY 


Event 

Date 

Time 

(GMT) 

Earthquake 

8-9-67 

4  3:26:50 

Zaza 

9-27-67 

17:00:00 

Earthquake 

10-4-67 

10:22:30 

Lanpher 

10-18-67 

14:30:00 

Rickey 

6-15-68 

16:00:00 

Sled 

8-29-68 

22:45:00 

Purse 

5-7-69 

13:45:00 

Hutch 

7-16-69 

14:55:00 

Pipkin 

10-8-69 

14:30:00 

Latitude 

CN) 

Longitude 

CW) 

Depth 

(km) 

Distance 
(degs 
to  LASA) 

39.9 

104.7 

5 

7 

37.1 

116.1 

0 

12 

38.5 

112.1 

18 

9 

37.1 

116.1 

0 

12 

37.3 

116.4 

0 

12 

37.3 

116.4 

0 

12 

37.3 

116.4 

0 

12 

37.1 

116.1 

0 

12 

37.3 

116.4 

0 

12 

Azimuth 

Location 

Magnitude, 

(NOS) 

170 

Denver 

5.3 

221 

Yucca 

- 

210 

Utah 

5.2 

221 

Yucca 

- 

223 

Pahute  Mesa 

- 

223 

Pahute  Mesa 

- 

223 

Pahute  Mesa 

5.8 

221 

Yucca 

5.6 

223 

Pahute  Mesa 

5.5 

TABLE  1-3 


DATA  BASE  FOR  NORSAR  RAYLEIGH-WAVE  STUDY 

Distance  Magnitude,  m, 

T  i  M  L  (  rf  n  A  rr  i  w\  -  _ _ 


Date 

Time 

(GMT) 

Latitude 

CN) 

Longitude 

CE) 

Depth 

(km) 

(degs 

to  NORSAR) 

Azimuth 

(deg) 

Location 

NOS 

u 

NORSAR 

6-15-71 

22:15:50 

41.4 

79.2 

N 

44 

85 

Kirgiz,  SSR 

5.1 

4.4 

7-11-71 

20:43:34 

37.2 

36.8 

9 

29 

131 

Turkey 

5.2 

4.4 

12-7-71 

12:04:57 

55.1 

-54.4 

N 

34 

291 

E.  of  Labrador 

5.4 

4.6 

12-16-71 

18:39:55 

77.9 

17.8 

N 

18 

5 

Svalbard  Region 

5.0 

4.2 

12-22-71 

07:05:24 

47.9 

48.2 

0 

25 

102 

Western  Kazakh 

6.0 

5.6 

3-21-72 

23:12:06 

35.8 

14.8 

N 

30 

171 

Mediterranean  Sea 

4.6 

4.4 

TABLE  1-4 

APPARENT  PHASE  VELOCITIES 


Period 

(sec) 


Phase  Velocity  (km/sec) 
LASA  NORSAR 


9 

3.11 

11 

3.25 

12 

3.26 

14 

3.29 

16 

3.27 

18 

3.35 

20 

3.56 

22 

3.57 

24 

3.62 

25 

- 

30 

3.70 

33 

- 

40 

3.85 

3.45 

3.50 

3.50 

3.60 

3.60 


3.80 

3.80 

4.10 


between  provinces  change  with  depth  and  different  frequency  Rayleigh  waves  are  sensitive  to 
different  depths.  Figure  1-19  shows  the  azimuthal  data  as  a  function  of  period  for  the  7  NTS  ex¬ 
plosions.  The  two  lines  are  the  theoretical  great  circle  azimuths  to  LASA  from  Pahute  Mesa 
and  Yucca  Flat.  The  waves  arrive  about  10  degrees  north  of  the  great  circle  path  at  short 
periods,  swinging  gently  to  the  south  until  the  azimuth  is  about  10  degrees  to  the  south  of  the 
theoretical  path  at  longer  periods.  Such  behavior  is  reasonable  for  a  path  that  crosses  a  transi¬ 
tion  zone  between  two  provinces,  where  the  crust  has  an  opposite  phase  velocity  contrast  for 
short  period  waves  compared  to  longer  period  waves.  In  Fig.  1-17  it  is  seen  that  the  mid-continent 
and  basin- range  phase  velocity  curves  merge  at  short  periods  with  the  mid-continent  values 
higher  at  longer  periods,  thus  the  longer  period  waves  should  be  pulled  to  the  south  as  observed. 

While  lateral  refraction  of  surface  waves  indicates  laterally  heterogeneous  structure  along 
the  propagation  path,  the  variation  of  phase  velocity  as  a  function  of  the  azimuth  of  an  event  indi¬ 
cates  heterogeneity  beneath  the  array.  Figure  1-20  plots  the  azimuth  of  the  arrival  of  the  20- sec 
component  of  the  Rayleigh  waves  for  the  NORSAR  events.  Figure  1-21  shows  the  measured  phase 
velocity  of  20-  and  16- sec  waves  for  the  various  azimuths.  The  peak  in  the  power  spectrum  of 
all  the  events  is  near  20  sec,  so  that  windowing  has  not  influenced  those  measurements.  The 
accuracy  of  the  data  is  indicated  by  the  scatter  of  the  values  at  16  sec.  Thus,  there  appears  to 

be  a  P-SV  velocity  anisotropy  beneath  the  array.  The  amount  is  about  5  to  10  percent,  which  is 
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comparable  to  the  7  percent  found  by  Wickens  and  Pec  and  Boore  for  40- sec  Love  waves  (SH) 
on  the  Canadian  shield.  The  most  likely  cause  of  anisotropy  is  the  heterogeneous  crustal  struc¬ 
ture  beneath  NORSAR.1  0 

A  final  observation  on  Rayleigh  waves  concerns  multipathing  of  the  events  recorded  at 
NORSAR.  Two  of  the  events  were  scanned  at  4  min  after  the  initial  Rayleigh  wave  arrival.  The 
wavenumber  spectra  are  given  in  Figs.  1-22  and  -23.  The  first  part  of  each  figure  shows  the 
initial  arrival.  The  latter  parts  show  later  arrivals  coming  in  at  various  azimuths.  The  arrival 
coming  from  the  north  in  Fig.  1-22  cannot  be  related  to  any  crustal  discontinuity.  The  LASA  and 
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NORSAR  bulletins  contain  no  event  which  might  have  produced  the  wave.  Therefore,  it  probably 
comes  from  the  original  event.  Its  genesis  is  puzzling.  The  other  arrival  can  be  related  to  the 
continental  margin  off  Britain,  and  the  later  arrival  in  Fig.  1-23  can  be  associated  with  the 
Black  Sea  depression.  A  note  of  interest  is  that  no  multipaths  were  observed  for  periods  longer 
than  18  sec,  implying  that  the  phenomenon  is  associated  with  upper  crustal  heterogeneity. 

T.  Landers 


F.  AN  EQUIVALENT  SCATTERING  SURFACE 

It  is  of  considerable  importance  to  explain  the  structure  within  the  earth  which  is  causing 

\ 

the  P-wave  amplitude  and  slowness  anomalies  presented  in  the  last  SATS.  One  of  the  simplest 
explanations  is  to  attribute  the  cause  of  the  anomalies  to  a  single  complex,  or  corrugated,  inter¬ 
face,  which  is  postulated  to  exist  in  the  vicinity  close  to  LASA.  Such  an  approach  has  been 
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taken  in  a  simple  form  by  Greenfield  and  Sheppard,  as  well  as  by  several  other  workers. 

Another  approach  entails  the  description  of  the  crust  and  upper  mantle  structure  under  LASA  as 
a  random  medium  which  can  be  characterized  by  a  correlation  distance,  as  discussed  in  detail  in 
Sec.  II- B.  Nevertheless,  there  is  a  certain  simplicity  involved  in  the  description  by  means 
of  a  single  corrugated  interface,  even  though  the  random  medium  approach  appears  preferable, 
if  we  try  to  account  both  for  amplitude  and  phase  anomalies.  In  addition,  there  are  several  very 
interesting  properties  that  can  be  ascribed  to  such  an  equivalent  scattering  surface.  It  is  for 
these  reasons  that  a  discussion  of  an  equivalent  scattering  surface  is  now  given.  This  discussion 

also  helps  to  link  the  present  work  with  that  which  has  been  done  previously  by  other  workers, 
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such  as  Greenfield  and  Sheppard. 

It  was  observed  that  the  average  curves  for  azimuth  deviation  (AZDEV),  and  the  slowness 
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correction  (SLOCOR)  data  in  Figs.  Ill- 3(a)  and  (b),  respectively,  resemble  the  results  that  would 
be  obtained  for  a  single,  or  possibly  multiple-parallel,  plane-dipping  interface.  Thus,  an  at¬ 
tempt  was  made  to  find  a  plane-dipping  interface  which  would  produce  values  for  AZDEV  and 
SLOCOR  that  would  fit  best  the  average  values  found  in  this  figure.  This  fit  was  performed  by 
using  a  least- squares  procedure  as  follows.  The  dip  angles  (DIP)  and  the  azimuths  of  the  up- 
dip  direction  (AZUPDIP)  for  a  set  of  plane-dipping  interfaces  were  specified,  as  well  as  the 
velocity  ratio,  which  was  allowed  to  vary  between  0.70  and  0.95;  i.e.,  V^/V^  varied  between 
these  limits,  where  V^  and  V^  are  the  compressional  velocities  in  the  crust  and  mantle,  re¬ 
spectively.  For  each  DIP  and  AZUPDIP,  and  for  each  event,  for  which  the  true  azimuth  and 
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angle  of  incidence  in  the  mantle  are  known,  the  computational  procedure  of  Niazi  was  used  to 
find  AZDEV  and  SLOCOR.  These  values  were  then  subtracted  from  the  corresponding  average 
values,  shown  in  the  aforementioned  figure,  for  each  event.  A  sum  of  squares  of  such  differ¬ 
ences  was  found  by  summing  the  squares  of  these  differences  for  each  of  the  31  events,  with  the 
SLOCOR  values  scaled  up  by  a  constant  multiplier,  so  that  they  are  comparable  to  the  AZDEV 
values.  The  best  fitting  plane  interface  is  defined  as  the  one  with  the  DIP  and  AZUPDIP  for 
which  the  minimum  sum  of  squares  is  obtained.  The  results  of  this  computation  are  shown  in 
Fig.  I-24(a)  for  a  velocity  ratio  of  0.75.  In  this  figure,  the  contours  of  constant  sum  of  squares 
are  shown  vs  DIP  and  AZUPDIP.  The  level  of  each  contour  is  given  in  dB  relative  to  the  min¬ 
imum  value  of  the  sum  of  squares.  Thus,  the  point  labeled  zero  dB  in  Fig.  I- 24(a)  corresponds 
to  the  best  fitting  plane  interface,  for  which  we  find  DIP  =  11  degrees  and  AZUPDIP  =  150  degrees. 
The  manner  in  which  AZDEV  and  SLOCOR  vary  as  a  function  of  azimuth  for  the  31  events,  for 
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the  best  fitting  plane  interface,  is  shown  in  Figs.  I- 24(b)  and  (c),  respectively.  It  is  interesting 
to  compare  these  results  with  those  in  Figs.  Ill- 3(a)  and  (b),  respectively,  in  Ref.  1. 

Another  interesting  observation  is  that  the  DIP  and  AZUPDIP  for  the  best  fitting  plane  inter¬ 
face  agree  very  well  with  the  corresponding  parameters  for  the  corrugated  surface  proposed  by 
Greenfield  and  Sheppard  for  a  velocity  ratio  of  0.75.  In  fact,  the  direction  of  strike  for  both 
sets  of  results  is  the  same.  This  agreement  is  quite  remarkable,  and  encouraging,  in  view  of 
the  different  measurement  techniques  which  were  employed  in  each  case. 

A  best  fitting  plane  interface  was  found  for  other  velocity  ratios  and  the  results  are  given 
in  Table  1-5.  The  AZUPDIP  remains  approximately  the  same  for  all  velocity  ratios  considered, 
while  the  DIP  increases  rapidly  as  the  velocity  ratio  approaches  unity.  This  latter  result  is  to 

be  expected.  There  is,  of  course,  no  way  that  the  velocity  ratio  can  be  known  with  any  certainty. 

13 

However,  some  seismic  refraction  work  done  by  Borchert  and  Roller  at  LASA  indicates  that 
a  velocity  contrast  of  0.75  is  a  reasonable  assumption. 

TABLE  1-5 

PARAMETERS  OF  BEST  FITTING,  PLANE- DIPPING 
INTERFACE  FOR  VARIOUS  VELOCITY  RATIOS 


Velocity 

Ratio 

Dip  Angle 
(deg) 

Azimuth  of  Up- Dip 
Direction  (deg) 

0.70 

9 

150 

0.75 

11 

150 

0.80 

13 

160 

0.85 

19 

160 

0.90 

27 

160 

0.95 

44 

160 

The  results  presented  previously  for  AZDEV  and  SLOCOR  were  used  to  compute  the  dip 
angles  and  azimuths  of  the  up-dip  direction,  for  a  set  of  plane-dipping  interfaces,  each  of  which 

is  assumed  to  be  responsible  for  the  measured  anomalies  at  each  subarray  and  for  each  event. 

1  2 

The  computational  procedure  employed  is  new  and  represents  an  inverse  to  Niazi's  method. 

The  computations  were  performed  for  various  velocity  ratios  between  0.70  and  0.95;  i.  e., 

VC/VM  varied  between  these  limits.  It  should  be  noted  that  it  was  not  always  possible  to  find 
a  solution  for  DIP  and  AZUPDIP  for  each  pair  of  anomalies  (AZDEV  and  SLOCOR)  which  were 
measured.  However,  it  was  possible  to  find  such  a  solution  in  the  great  majority  of  cases. 

The  DIP  and  AZUPDIP  obtained  for  a  velocity  ratio  of  0.75  are  shown  plotted  vs  the  true 
azimuth  in  Figs.  I- 25(a)  and  (b),  respectively.  In  addition,  a  mean,  or  average,  curve  is  drawn 
in  these  figures,  where  the  average  is  taken  over  the  nine  values  obtained  within  the  subarrays 
for  each  event.  The  average 'curves  for  DIP  and  AZUPDIP  in  Figs.  I-25(a)  and  (b)  vary  quite 
randomly  about  an  overall  average  of  25.1  and  157.4  degrees,  respectively.  Some  values  of  DIP 
can  be  as  high  as  75  degrees.  In  addition,  there  does  not  seem  to  be  any  correlation  between 
the  values  of  DIP  and  AZUPDIP. 

It  has  been  shown  that  the  values  for  DIP  and  AZUPDIP,  corresponding  to  the  AZDEV  and 
SLOCOR  data,  tend  to  scatter  considerably  about  their  mean  values.  This  makes  it  difficult  to 
discuss  these  quantities  in  a  deterministic  sense.  It  appears  more  useful  to  treat  DIP  and 
AZUPDIP  as  random  variables  which  assume  values  in  an  appropriate  sample  space.  The  sample 
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spaces  for  DIP  and  AZUPDIP  consist,  of  course,  of  the  real  line  between  0  and  90,  and  the  real 

line  between  0  and  360  degrees,  respectively.  We  now  show  how  to  compute  the  theoretical 

probability  distribution  functions  (PDF)  for  DIP  and  AZUPDIP,  and  how  these  compare  with  the 

empirical  distribution  functions  (EDF)  for  these  quantities. 

In  order  to  derive  the  PDFs  for  DIP  and  AZUPDIP,  we  will  use  some  well-known  results 

1 4 

developed  by  Longuet- Higgins  in  the  study  of  random  sea  waves,  or  random  surfaces,  as  well 

1 5 

as  results  due  to  Middleton.  It  was  assumed  that  the  particular  structure  under  LASA  that  is 
causing  the  P-wave  anomalies  is  a  complicated  corrugated  surface  which  can  be  described  as 
the  sum  of  a  planar  surface  plus  a  random  surface.  This  planar  surface  corresponds  to  the  best 
fitting  plane  interface  described  previously.  The  random  surface  is  assumed  to  be  a  sample 

1 4 

function  from  an  isotropic,  Gaussian,  random  process  with  zero  mean,  cf. ,  Longuet- Higgins. 

1 4 

Using  these  assumptions,  as  well  as  some  results  of  Longuet- Higgins  and  Middleton  (pp.  411- 
1 5 

418)  it  was  shown  that  the  probability  density  functions  for  the  local  slope  a  and  AZUPDIP  9 
of  the  random  surface  are  given  by 

pg(  a )  =  (2<*/M)  exp[-(a2  +  a2)/M]  ,  a  >0 

=  0  ,  a  <  0 

pA(0  )  =  (l/27r)  exp  (-  a  2)  ({l  +  nTJt  a  cos  (9-9  )  exp  [a  2  cos2(9  -  9  )]} 

X  { 1  +  erf  [aQ  cos  (0  -  0Q)]}^  ,  0  ^  9  ^  2tt 

where  IQ(x)  is  the  zero-order  modified  Bessel  function  of  the  first  kind,  erf(x)  is  the  error  func- 
ti°n,_ja0  and  eo  are  the  sloPe  and  AZUPDIP,  respectively,  for  the  best  fitting  plane  interface, 

M  =  or  —  and  aQ  =  or  /M  where  a  denotes  the  expectation,  or  average,  of  the  random  vari¬ 
able  o?2.  The  PDFs  for  a  and  9  can  be  obtained  from  the  probability  density  functions  given 

above  in  the  usual  manner.  It  should  also  be  noted  that  if  a  =0,  the  PDFs  for  a  and  9  reduce 

°  16 

to  the  Rayleigh  and  uniform  distributions,  respectively,  cf. ,  Davenport  and  Root.  The  local 

_  i 

dip  6  is  related  to  the  local  slope  a  in  the  standard  way,  i.  e. ,  6  -  tan  (o?).  Thus,  the  PDF 
for  6  can  also  be  found  from  the  PDF  for  a,  but  it  is  convenient  to  work  only  with  the  PDF 
for  a. 

The  EDFs  for  a  and  9  were  computed  from  the  AZDEV  and  SLOCOR  data  using  the  proce¬ 
dure  mentioned  previously,  and  the  results  are  shown  in  Figs.  I- 26(a)  and  (b),  respectively,  for 

2 

a  velocity  ratio  of  0.75.  The  EDF  for  a  can  be  used  to  estimate  a  ,  i.e.,  its  mean-square 
value,  from  which  the  value  for  M  and  aQ  can  be  obtained  using  the  above  equations,  since 
is  known,  i.e.,  oiQ  =  tan  (11°  )  =  0.19,  cf..  Table  1-5.  The  PDFs  for  a  and  0  can  then  be  com¬ 
puted  using  the  equations  given  previously,  and  the  results  are  given  in  Figs.  I- 26(a)  and  (b).  It 
is  seen  from  this  figure  that  there  is  reasonably  good  agreement  between  the  EDFs  and  PDFs. 

This  means  that  the  DIP  and  AZUPDIP  results  obtained  from  the  AZDEV  and  SLOCOR  data  at 
LASA  correspond  closely  to  the  results  which  would  be  obtained  if  one  were  to  measure  the  local 
DIP  and  AZUPDIP  of  a  Gaussian  random  surface  which  is  superimposed  on  a  planar  surface. 

An  attempt  was  made  to  describe  graphically  the  random  surface  discussed  previously,  and 
the  computations  employed  are  described  briefly.  A  velocity  ratio  of  0.75  was  chosen  and  a 
60-  x  60-km  region  placed  symmetrically  about  the  A0  subarray  at  LASA  was  considered.  A 
master  plane  was  taken  passing  arbitrarily  at  50  km  beneath  subarray  A0  with  a  DIP  and  AZUPDIP 
of  11  and  150  degrees,  respectively.  This  plane  corresponds  to  the  best  fitting  plane-dipping 
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interface  discussed  previously.  The  direction  cosines  of  the  seismic  rays  in  the  crust  were 
computed  from  the  AZDEV  and  SLOCOR  data  using  the  procedure  mentioned  previously.  These 
rays  were  then  traced  back  to  their  point  of  intersection  with  the  master  plane.  In  the  vicinity 
of  these  points,  the  master  plane  was  modified  to  have  a  local  DIP  and  AZUPDIP,  which  are 
computed  from  the  AZDEV  and  SLOCOR  data.  The  resulting  surface  was  then  smoothed  once 
with  a  standard,  five-point  smoothing  scheme.  This  smoothed  surface  is  shown  in  Fig.  1-27, 
which  gives  the  contours  of  constant  depth  (in  km)  as  a  function  of  the  x-  and  y- coordinates.  A 
three-dimensional  view  of  the  surface  is  depicted  in  Figs.  I-28(a-d).  The  viewing  angle  is 
30  degrees,  relative  to  the  vertical  z-axis,  and  the  rotation  angles  are  30,  120,  210  and  300 
degrees,  as  indicated  in  Figs.  I-28(a-d),  respectively. 

J.  Capon 
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Fig.  1-1.  P-wave  amplitudes  at  NORSAR  subarrays  for  a  presumed  explosion 
in  Eastern  Kazakh,  25  May  1971.  SNR  at  2  Hz  is  also  shown. 


217m/i./sec 


4C  BEAM 


lH-2-10T7?| 


mu  /car 


y\[\f\J\ 


—  -  - -I —  - ^ — 'N-'‘ 


Fig.  1-2.  P- waves  at  sensors  of  two  high  amplitude  NORSAR  subarrays  for  presumed 
explosion  of  25  May  1971  in  Eastern  Kazakh.  Subarray  beams  are  also  displayed. 
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Fig.  1-3.  P- waves  at  sensors  of  two  low- amplitude  NORSAR  subarrays 
for  same  event  as  in  Fig.  1-2. 
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Fig.  1-4.  NORSAR  steered  beams  of  seven 
events  using  center  sensors.  Events  3  and  4 
are  presumed  explosions  in  Eastern  Kazakh. 
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Fig.  1-6.  NORSAR  amplitude  for  event  2.  Azimuth  =  71.8°,  A  =  67.1°. 
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Fig.  1-7.  NORSAR  amplitude  for  event  3.  Azimuth  =  73.2°,  A  =  39. 


Fig.  1-8.  NORSAR  amplitude  for  event  4.  Azimuth  =  74.0°,  A  =  38. 
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Fig.  I- 10.  NORSAR  amplitude  for  event  6.  Azimuth  =  93.7°,  A  =  44.8°. 
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Fig.  I- 1 1 .  NORSAR  amplitude  for  event  7.  Azimuth  =  94.7°,  A  =  44.9°. 


Fig.  1-12.  dT/dA  values  measured  at  NORSAR  (crosses)  together  with  data 
of  Johnson2  and  calculated  curve  for  model  WNA  fitted  to  Johnson’s  data  and 
nuclear  explosion  travel  times  within  western  North  America. 
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Fig.  1-13.  Beams  formed  using  a  single 
NORSAR  subarray  (2B)  of  the  P  wave 
from  a  deep  Hindu  Kush  earthquake:  (a) 
with  hand-picked  time  delays;  (b)  with 
delays  calculated  for  a  plane  wave. 
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Fig.  1-14.  Distribution  in  slowness  space  of  events  used  to  measure  sensor  delays. 
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NUMBER  OF  SENSORS 


(a) 


Fig.  1-15.  Expected  spectral  degradation 
as  a  function  of  rms  time  error  for  per¬ 
fectly  coherent  signals.  (E  means  expec¬ 
tation;  f  is  frequency;  AT  is  the  timing 
error;  P(f)  is  the  "transfer  function"  for 
beamforming. ) 


(b) 


Fig.  1-16.  Distribution  of  rms  time  errors,  (a)  Observed  distribution 
for  individual  NORSAR  sensors,  (b)  Distribution  to  be  expected  if  a  con¬ 
stant  time  correction  is  used. 
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Fig.  1-17.  Phase  velocity  curve  for  LASA. 
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Fig.  1-18.  Phase  velocity  curve  for  NORSAR. 
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Fig.  1-19.  Azimuth  of  arrival  of  Rayleigh  waves  at  LASA  from  NTS 
as  a  function  of  period. 
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Fig.  1-20.  Azimuth  of  arrival  of  20-sec 
component  of  Rayleigh  wave  for  each 
NORSAR  event. 
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Fig.  1-21.  Phase  velocity  for  16-  and  20-sec  Rayleigh  waves  at  NORSAR 
as  a  function  of  azimuth. 
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Fig.  1-22.  Wavenumber  spectra  for  14-  and  48-sec  Rayleigh  waves  at  NORSAR 
from  event  of  7  December  4971.  Plot  (a)  shows  main  arrival;  plots  (b)  and  (c), 
for  data  4  minutes  later,  show  multipath  arrivals. 
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Fig.  1-23.  Wavenumber  spectra  for  event  of  11  July  1971  at  NORSAR. 
Plot  (a)  shows  main  arrival;  plot  (b),  for  data  3  min.  10  sec  later, 
shows  a  multipath. 
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VELOCITY  RATIO  «  0.75 
DIP  ANGLE  ■  IT 

AZIMUTH  OF  UP-DIP  DIRECTION  «  150* 


Fig.  1-24. 


Determination  of  parameters  of  best  fitting  plane  dipping  interface. 
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Fig.  1-25.  Dip  angle  and  azimuth  of  up-dip  direction  vs  true  azimuth. 
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AZIMUTH  OF  UP-DIP  DIRECTION  »  150* 


Fig.  1-26.  Comparison  of  empirical  distribution  functions  of  slope  and  azimuth 
of  up-dip  direction  with  theoretical  cumulative  distribution  functions. 
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Velocity  Rotio  ■  0.75 
Dip  Angie  *  11° 

Azimuth  of  Up-Dlp  Direction  »  150° 


Fig.  1-27.  Contour  plot  of  depth  vs  x-  and  y-coordinates  for  equivalent 
scattering  surface. 
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VELOCITY  RATIO  •  0.75 
VIEW  ANGLE  -  30*  FROM  VERTICAL 


Fig.  1-28.  Three-dimensional  views  of  equivalent  scattering  surface. 
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II.  EARTH  STRUCTURE 


A.  REPRESENTATION  OF  APPROACH  ANGLE  ANOMALIES 

Approach  angle  anomalies  have  been  measured  for  P-wave  signals  at  LASA  observed  during 
the  first  5  sec  after  the  onset  time.  These  anomalies  were  obtained  within  the  nine  subarrays 
comprising  the  A-,  B-,  and  C-rings  at  LASA,  and  the  results  were  presented  in  the  previous 

A 

SATS.  Azimuth  deviation  (AZDEV)  and  slowness  (or  dT/dA)  correction  (SLOCOR)  data  were 

plotted  vs  the  true  azimuth  of  the  epicenter  of  the  event.  It  was  pointed  out  that  this  type  of 

presentation  of  the  results  held  certain  advantages  in  the  interpretation  of  the  data. 

We  can,  however,  represent  the  anomalies  in  a  circular  diagram  in  a  way  used  by  Manchee 
2  3 

and  Weichert  and  Vinnik  and  Nikolayev.  This  approach  entails  the  plotting  of  the  anomalies  as 

a  series  of  arrows  in  (dT/dA,  Azimuth)  space.  The  tail  of  the  arrow  represents  the  dT/dA  and 

azimuth  measured  for  the  incoming  P-wave  signal  at  the  array,  and  its  head  represents  the 

dT/dA  and  azimuth  that  would  be  obtained  from  the  true,  or  in  this  case  CGS,  location  and  the 

4 

J-B  tables.  This  method  has  been  applied  recently  by  Davies  and  Sheppard  to  the  study  of  the 
P-wave  slowness  anomalies  at  LASA.  They  have  termed  the  resultant  presentation  the  "array 
diagram." 

2 

It  was  shown  by  Manchee  and  Weichert,  that  if  the  anomaly  were  caused  by  a  plane-dipping 
interface  located  directly  under  the  array,  the  array  diagram  would  consist  of  a  set  of  error 
vectors  whose  magnitudes  are  all  approximately  the  same  and  which  all  point  in  the  direction  of 
the  dip  of  the  interface.  Thus,  the  effects  of  a  sloping  interface  can  be  determined  from  the  array 
diagram  by  fitting  a  constant  vector  to  the  error  vectors  in  the  diagram. 
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Additional  results  in  the  interpretation  of  the  array  diagram  are  due  to  Davies  and  Sheppard. 

They  point  out  that  a  smoothly  varying  distribution  of  arrows,  all  approximately  of  equal  length 

and  pointing  in  the  same  direction,  would  indicate  a  bias  introduced  beneath  the  array.  The 

reason  for  this  is  that  all  ray  paths  from  the  various  event  epicenters  to  LASA  must  traverse 

the  same  structure  under  the  array.  This  result  is  slightly  different  from  that  of  Manchee  and 
2 

Weichert,  mentioned  previously,  in  that  the  anomaly  caused  by  a  smooth  distribution  of  arrows, 
although  due  to  structure  under  the  array,  need  not  be  caused  by  a  plane-dipping  interface.  In 
addition,  Davies  and  Sheppard  point  out  that  the  spatial  rate  of  variation  of  the  error  vectors  can 
be  associated  with  heterogeneities  at  certain  depths  in  the  earth.  The  more  rapid  these  spatial 
variations  become,  the  more  distant  from  LASA  the  heterogeneities  must  be  placed  to  explain 
the  anomalies.  The  reason  for  this  can  be  easily  seen  on  the  basis  of  the  separation  of  the  ray 
paths  from  the  various  event  epicenters  to  LASA. 

In  view  of  the  power  of  the  presentation  contained  in  the  array  diagram,  the  P-wave  anom¬ 
alies  measured  within  subarrays  at  LASA  have  been  displayed  in  "subarray  diagrams"  and  the 
results  for  the  B-ring  are  given  in  Figs.  Il-l(a-d).  The  four  diagrams  in  Fig.  II -1  are  typical 
of  the  nine  subarray  diagrams  obtained.  There  are  certain  regularities  seen  both  within  each 
diagram  and  amongst  the  several  diagrams.  For  instance  in  all  diagrams  there  is  a  tendency 
for  the  arrows  to  line  up  in  a  NW-SE  direction.  This  suggests  that  dipping  structures  of  some 
sort  underlie  each  subarray,  and  that  these  structures  are  moderately  persistent  across  the 
whole  of  the  C-ring.  This  is  shown  quite  clearly  for  instance  by  the  similarity  of  all  diagrams 
to  a  subarray  average  diagram  [Fig.  II-2 (a)  ],  being  the  average  of  nine  subarrays.  Figure  II- 2 (b) 
shows  the  mean  vector  of  Fig.  II-2(a)  applied  uniformly  to  all  sources,  and  Fig.  II-2(c)  shows  the 
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theoretical  effect  of  a  single  dipping  interface  with  a  strike  of  46°  to  226°  and  a  dip  of  1 1  °  (veloc¬ 
ity  contrast  at  interface  is  0.75:1).  This  is  used  only  as  an  example.  Niazi's  computational  pro- 
5 

cedure  has  been  used.  The  result  agrees  quite  favorably  with  the  corresponding  results  of  1 1  0 
and  330°  found  by  an  alternative  method,  cf.,  Sec.I-F. 

4 

Davies  and  Sheppard  found  that  the  mean  error  vector  for  their  array  diagram  pointed  north 
and  had  a  magnitude  of  0.25  sec/deg.  It  is  possible  to  show  that  this  error  vector  could  be  as¬ 
cribed  to  a  plane-dipping  interface  under  LASA  with  dip  of  about  3°  and  direction  of  dip  at  an 
azimuth  of  0°. 

It  is  encouraging  to  note  that  there  are  certain  similarities  between  the  subarray  average 
diagram  given  in  Fig.II-2(a)  and  the  array  diagram  for  the  full  LASA  given  in  Fig.  1  of  Ref.  4. 
However,  the  fact  that  there  are  differences  can  be  explained  by  the  following  reasons: 

(1)  The  effect  of  the  random  part  of  the  local  crust  and  upper  mantle  structure  at  LASA 
is  greater  in  the  subarray  average  diagram  than  in  the  LASA  diagram  (cf..  Sec.  II-C). 

(2)  The  subarray  average  diagram  represents  an  average  taken  over  a  circular  area 
of  the  Moho  of  about  30  km  in  diameter,  while  the  array  diagram  averages  over  a 
circle  of  200  km  diameter.  Thus,  the  effect  of  a  possibly  warped  Moho  (if  the  Moho 
is  the  prime  interface)  under  LASA  may  be  seen  in  the  subarray  average  diagram, 
but  not  in  the  array  diagram. 

(3)  The  array  diagram  is  based  on  a  least-squares  fit  of  a  plane  wave  to  arrival  times 

measured  from  the  first  motion  of  the  P-wave.  As  a  result,  this  diagram  pertains 
to  the  high-frequency  part  of  the  P-wave  spectrum.  The  subarray  average  diagram 
was  obtained  at  a  frequency  of  0.8  Hz.  In  addition,  there  may  be  some  effects  due 
to  frequency  windowing  in  the  results  reported  here  that  do  not  appear  in  the  array 
diagram.  j.  Capon 

D.  Davies 

B.  CHARACTERIZATION  OF  CRUST  AND  UPPER  MANTLE  STRUCTURE 

UNDER  LASA  AS  A  RANDOM  MEDIUM 

An  analysis  has  been  made  of  the  P-wave  slowness  and  amplitude  anomalies  within  the  sub- 

1 

arrays  at  LASA.  The  method  of  measurement  was  described  in  the  last  SATS.  These  anomalies 
were  measured  for  the  P-waves  of  31  events  located  at  various  azimuths  and  distances  from 
LASA.  They  were  obtained  within  the  subarrays  contained  in  the  A-,  B-,  and  C -rings  at  LASA, 
so  that  a  total  of  279  observations  were  taken.  The  data  were  observed  on  the  short-period  ver¬ 
tical  seismometers  at  LASA.  The  frequency  used  in  the  measurement  was  0.8  Hz,  so  that  the 
wavelength  of  the  P-wave  is  about  10  km,  assuming  the  average  compressional  velocity  is  about 
8  km/sec. 

A  model  has  been  proposed  to  explain  the  cause  of  the  P-wave  anomaly  data,  in  which  the 
crust  and  upper  mantle  structure  under  LASA  is  characterized  as  a  random,  or  Chernov-type, 
medium.^  This  model  seems  to  be  quite  successful  in  explaining  the  P-wave  anomaly  data  and 
a  description  of  these  results  is  now  presented.  In  order  to  do  this,  it  is  necessary  to  present 
a  brief  review  of  some  of  the  results  due  to  Chernov  concerning  wave  propagation  in  a  random 
medium.  Although  these  results  are  derived  for  the  case  of  wave  propagation  in  a  liquid,  or 
gaseous,  medium,  it  may  be  shown  that  they  can  be  used  for  an  elastic  medium  with  only  minor 
modifications. 
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The  index  of  refraction  for  the  random  medium  is  the  random  function  p(x ,  y,  z)  defined  as 

H(x,  y,  z)  =  [c(x,  y,  z)/cq(x,  y,  z)  ]  -  1 

where  c(x,  y,  z)  is  the  velocity  and  cq(x,  y,  z)  is  the  average  velocity,  i.e.,  cq(x,  y,  z)  =  c  (x,  y,  z). 

In  addition,  it  is  assumed  that  only  small  variations  from  homogeneity  exist;  i.e.,  the  mean- 
square  value  of  the  relative  index  of  refraction  is  much  less  than  unity 

(jl2(x,  y,  z)  «  1 

If  the  random  process  is  statistically  isotropic,  then  the  correlation  coefficient,  N,  depends  only 
on  the  distance,  r,  between  the  points  (x^,  y^ ,  z^),  (x^,  y^,  z^).  A  particularly  useful  correlation 
coefficient,  which  was  assumed  for  the  random  medium  in  the  analysis,  is  the  Gaussian  correla¬ 
tion  coefficient,  given  by 

N(r)  =  exp(-r2/a2)  , 

where  a  is  known  as  the  correlation  distance. 

The  propagation  of  a  wave  in  a  medium  with  random  inhomogeneities  leads  to  fluctuations 
of  the  amplitude  and  phase  characteristics  of  the  wave  field  due  to  the  superposition  of  the  scat¬ 
tered  waves  and  the  primary  wave.  There  is  a  dependence  between  the  fluctuations  of  the  char¬ 
acteristics  of  the  wave  field  and  the  fluctuations  of  the  relative  index  of  refraction.  Thus,  it  is 
feasible  to  relate  the  statistical  properties  of  the  wave  field  to  the  statistical  properties  of  the 
ensemble  of  random  media.  Hence,  it  is  possible,  by  measuring  the  amplitude  and  phase  char¬ 
acteristics,  to  determine  several  parameters  which  are  important  in  the  description  of  the  ran¬ 
dom  medium.  A  comprehensive  set  of  such  results  has  been  given  by  Chernov.^  The  particular 
model  assumed  by  Chernov  is  shown  in  Fig.  II- 3 .  A  uniform  plane  wave  is  propagating  in  the 
left  half-space,  x  <  0,  which  is  homogeneous,  i.e.,  contains  no  random  inhomogeneities.  This 
wave  is  termed  the  direct  or  homogeneous  wave,  and  enters  the  right  half-space,  x  >  0,  which 
is  inhomogeneous  and  contains  random  inhomogeneities.  The  wavenumber,  in  rad/km,  and  fre¬ 
quency,  in  rad/sec,  of  the  wave  in  the  left  half-space  are  denoted  by  k  and  co,  respectively. 

Thus,  we  have  k  =  cj/cq  =  2.n/\,  where  \  is  the  wavelength  in  km. 

It  is  assumed  that  there  are  two  receivers  located  at  the  coordinates  (L,  0,  0)  and  (L,  0,  l ) 
as  shown  in  Fig.  H-3.  An  essential  role  in  the  theory  is  played  by  the  wave  parameter  D  defined 
as:  D  =  4L/ka  .  It  is  assumed  that  large-scale  inhomogeneities  exist  in  the  random  medium, 
i.e.,  ka  »  1,  and  also  that  L  »  a. 

The  amplitude  fluctuation,  B,  is  defined  as:  B  =  ln(A/AQ),  where  Aq  is  the  amplitude  of  the 
wave  in  the  left  half- space  and  A  =  A(x,  y,  z)  is  the  amplitude  of  the  wave  field  at  the  point  (x,  y,  z) 
in  the  right  half-space.  The  phase  fluctuation,  S,  is  defined  as  the  phase  of  the  wave  field  in  the 
right  half-space  relative  to  that  of  the  homogeneous  wave.  The  ratio  of  the  mean-square  fluctua¬ 
tions  of  phase  and  amplitude,  B2/S2,  was  of  importance  in  the  analysis.  Using  the  results  of 
Chernov  (pp.  82-83),  we  have 

1  -  -{j  tan"1  (D) 

"^2  1  +  tan  1  (D) 

In  addition,  the  crosscorrelation  coefficient  between  the  amplitude  and  phase  fluctuations  was 
derived  by  Chernov  (pp.  85-89),  as 
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Rbg  =  BS/(B2S2)1/2 


=  l/2in(l  +  D2)/{D2  -  [tan  1  (D)]2}1//2 


The  autocorrelation  of  the  amplitude,  as  well  as  phase,  fluctuations  in  the  direction  which 
is  longitudinal  to  the  direction  of  wave  propagation  in  the  left  half-space,  i.e.,  in  the  direction 
of  the  x-axis,  was  also  computed  by  Chernov  (pp.  89-95).  These  correlations  were  computed 
between  two  receivers,  or  sensors,  located  at  the  coordinates  (L^,  0,  0)  and  (L^,  0,  0).  It  was 
shown  that  both  the  amplitude  and  phase  fluctuations  are  practically  completely  correlated  over 
a  longitudinal  distance,  AL,  given  by:  AL  =  ka  /2.  Since  it  has  been  assumed  that  ka  »  1,  it 
follows  that  AL  »  a. 

The  transverse  autocorrelation  coefficient  for  the  amplitude  and  phase  fluctuations,  com¬ 
puted  between  the  two  receivers  shown  in  Fig.  II- 3,  which  are  located  at  the  coordinates  (L,  0,  0) 
and  (L,  0,  i)  was  derived  by  Chernov  (pp.  95-110).  These  two  receivers  are  located  in  a  plane 
which  is  transverse  to  the  direction  of  propagation  of  the  wave  in  the  left  half- space,  i.e.,  the 
direction  of  the  x-axis.  The  results  for  the  case  when  the  wave  parameter  is  large,  D  »  1, 
i.e.,  Fraunhofer  diffraction,  are  given  by 


iyi)  =  b^/b2 


Rg(i)  =  S1S2/S2 


- 1 - ~\ - 

1  -  ^  tan  (D) 

1  +  ^  tan-1  (D) 


where  B^,  B^,  S^,  are  the  amplitude  and  phase  fluctuations  observed  at  the  two  receivers  and 
Si(-)  is  the  sine  integral  function  defined  as 

0.,  ,  f X  sin  (u) 

Si(x)  =  \  du 

o  u 


It  was  also  shown  by  Chernov  that  the  values  for  these  coefficients  are  near  unity  when  I  «  a, 
and  decrease  to  zero  when  l  »  a.  Although  this  result  is  derived  for  the  case  of  the  Gaussian 
correlation  coefficients,  it  has  been  shown  to  be  true  for  a  wide  class  of  correlation  coefficiency 
by  Tatarski  (p.  146).  Thus,  the  measurement  of  these  coefficients  can  provide  a  powerful 
method  for  determining  the  correlation  distance  of  the  random  medium,  which  is  an  important 
parameter. 

An  indication  is  now  given  of  the  manner  in  which  Chernov’s  theory  of  wave  propagation  in 
a  random  medium  was  used  to  characterize  random  properties  of  the  crust  and  upper  mantle 
under  LASA,  using  the  P-wave  anomaly  data  presented  previously.  Toward  this  end,  it  is  nec¬ 
essary  to  investigate  the  way  in  which  the  P-wave,  or  compressional,  velocity  varies  with  depth 
beneath  LASA.  A  series  of  seismic  refraction  profiles  have  been  made  at  LASA  by  Borcherdt 

and  Roller/*  Warren,  et  aL/*  and  Gibbs.10  In  addition,  there  are  studies  of  the  velocity  structure 

1 1 

in  the  upper  mantle  for  the  western  continental  United  States  made  by  Niazi  and  Anderson, 

12  13 

Archambeau,  et  aL,  and  Bolt  and  Nuttli.  If  an  attempt  is  made  to  combine  all  the  preceding 
results  about  the  velocity  structure  in  the  crust  and  upper  mantle,  a  velocity  structure  such  as 
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that  in  Fig.  II-4(a)  may  be  obtained.  The  hyperfine  variations  in  the  velocity  structure  have  been 
placed  there  in  order  to  explain  the  mechanism  which  causes  the  P-wave  anomaly  data  at  LASA, 
and  do  not  exist  in  any  of  the  preceding  works.  The  results  shown  in  Fig.  ll-4(a)  pertain  to  a 
particular  point  on  the  surface  of  the  earth,  and  the  velocity  structure  (as  a  function  of  depth) 
will  change  as  this  point  is  allowed  to  vary.  This  variation  between  two  sets  of  results  will  be 
small  if  the  points  are  within  a  correlation  distance  of  each  other.  This  correlation  distance 
was  measured  and  found  to  be  about  12  km,  as  will  be  discussed  subsequently. 

The  velocity  structure  in  Fig.  11-4  shows  the  velocity  varying  randomly  about  8  km/sec  in 
the  upper  mantle,  at  approximately  60-  to  200 -km  depth.  The  influence  of  the  low- velocity  zone 
at  depths  between  90  and  150  km  is  also  noted.  The  velocity  change  due  to  the  Mohorovicic  dis¬ 
continuity  probably  occurs  at  about  50  km,  while  there  is  another  sudden  change  at  about  20  km. 

In  addition,  there  is  a  large  change  in  velocity  that  occurs  about  5  km  beneath  the  surface  due 
to  sedimentary  layers. 

Unfortunately,  the  velocity  structure  shown  in  Fig.  II-4(a)  does  not  correspond  to  that  which 
would  be  expected  in  a  random  medium  which  is  spatially  homogeneous  and  statistically  isotropic. 
That  is,  the  correlation  function  for  the  index  of  refraction  of  such  a  medium  depends  on  position 
in  the  medium  and  not  just  on  differences  in  position,  as  discussed  previously.  There  appears 
to  be  no  simple  way  to  modify  Chernov’s  results  to  take  this  nonhomogeneous  and  non-isotropic 
condition  into  account.  Thus,  in  order  to  use  the  results  given  by  Chernov,  it  is  necessary  to 
make  some  simplifying  assumptions  about  the  velocity  structure,  but  yet  retain  the  essence  of 
this  variation.  Such  a  simplified  velocity  structure  is  shown  in  Fig.  II - 4(b) .  It  is  possible  to 
replace  the  velocity  structure  in  Fig.ll-4(b)  with  that  in  Figs.  II-4(c)  and  (d),  which  depict  the 
velocity  structure  for  a  random  and  a  deterministic  medium,  respectively.  In  other  words,  as 
a  gross  approximation,  it  is  possible  to  consider  the  amplitude  and  phase  fluctuations  produced 
by  propagation  through  the  medium  with  velocity  structure  given  in  Fig.  II-4(b)  as  equivalent  to 
those  produced  by  propagation  through  the  medium  whose  velocity  structure  is  given  in  Fig.  II-4(c), 
followed  by  the  medium  whose  structure  is  shown  in  Fig.  II-4(d).  An  assumption  has  also  been 
made  that  the  sedimentary  layers  play  a  relatively  minor  role  in  causing  the  P-wave  anomalies 
at  LASA.  This  is  a  simplifying  assumption  that  cannot  be  completely  justified  at  this  time. 

Thus,  the  problem  of  considering  the  effect  of  the  crust  and  upper  mantle  on  P-waves  was 
reduced  to  that  of  considering  the  random,  or  Chernov-type,  medium  in  Fig.ll-4(c),  followed 
by  the  deterministic  medium  in  Fig.  II-4(d).  The  effect  of  the  deterministic  medium  can  be  taken 
into  account  as  follows.  The  x-axis  in  Fig.  II- 3  is  oriented  so  that  it  lies  along  the  direction  of 
propagation  of  the  P-wave  defined  by  the  refraction  which  takes  place  in  the  deterministic  medium 
under  LASA,  as  well  as  any  other  deterministic  medium  external  to  the  LASA  which  lies  along 
the  propagation  path,  such  as  in  the  lower  mantle  (or  source)  region. 

In  order  to  use  the  results  of  Chernov  on  transverse  autocorrelation  coefficients  to  determine 
the  correlation  distance,  it  is  necessary  to  measure  the  amplitude  and  phase  fluctuations  at  the 
receivers  located  at  the  points  PI’  and  P2  depicted  in  Fig.  II- 3;  i.e.,  the  two  receivers  should  lie 
in  the  plane  which  is  transverse  to  the  direction  of  propagation  of  the  wave  in  the  left  half-space, 
or  x-axis.  However,  it  is  only  possible  to  have  receivers  located  at  the  surface  of  the  earth  at 
points  PI  and  P2,  where  PI  is  the  point  of  intersection  of  the  surface  with  the  line  parallel  to, 
and  i  km  above,  the  x-axis,  as  shown  in  Fig.  II-3.  It  can  be  shown  that  the  correlation  of  the 
fluctuations  between  PI  and  P2  is  a  good  approximation  for  that  between  PI'  and  P2,  since  the 
correlation  of  the  fluctuations  between  PI’  and  PI  is  close  to  unity,  as  indicated  previously. 
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It  is  now  necessary  to  show  how  the  amplitude  and  phase  fluctuations,  required  in  Chernov’s 
theory,  may  be  obtained  from  the  P-wave  anomaly  data.  In  addition,  it  is  also  necessary  to  show 
how  these  quantities  may  be  used  to  estimate  the  important  parameters  that  characterize  the  ran¬ 
dom  medium,  which  were  discussed  previously.  A  basic  assumption  is  made  that  the  P-waves 
for  each  event  enter  a  random  medium,  which  can  be  considered  to  be  an  independent  realization 
from  the  ensemble  of  random  media.  In  reality  there  is,  of  course,  only  a  single  random  medium 
under  LASA.  However,  this  assumption  is  valid  if  the  epicenters  of  the  events  are  sufficiently 
separated  from  each  other  so  that  the  P-wave  trajectories  for  any  given  event,  in  the  upper  man¬ 
tle,  intersect  a  fixed  horizontal  plane,  for  example,  about  200  km  beneath  LASA,  at  points  which 
are  more  than  a  few  correlation  distances  from  those  of  any  other  event.  Thus,  the  ensemble 
(or  statistical)  averages  may  be  performed  by  averaging  the  results  obtained  for  the  individual 
events. 

If  the  measurement  for  the  relative  power  at  the  i-th  subarray  and  j-th  event  is  denoted  by 
P..,  then  the  amplitude  fluctuation  is  defined  as 


B..  =  -c(P..  -  P.)f  j  =  l . 31  , 


i  =  l . 9  , 


where  c  =  ln(10)/20,  and  the  average  amplitude  fluctuation  is 


v(J  -,)/* 


This  definition  for  the  amplitude  fluctuation  was  found  to  give  the  best  agreement  of  the  experi¬ 
mental  results  with  the  theory  given  by  Chernov.  The  estimate  of  the  transverse  autocorrelation 

1  4 

coefficient  for  the  amplitude  fluctuations  is  obtained  in  the  usual  manner,  cf.,  Cramer. 

In  order  to  compute  the  phase  fluctuation,  it  is  necessary  to  know  the  phase,  or  time  of 
arrival,  of  the  diffracted  wave  at  the  receiver  relative  to  that  of  a  hypothetical  homogeneous 
wave.  The  arrival  time  of  this  homogeneous  wave  is  the  time  at  which  the  wave  in  the  left  half¬ 
space,  shown  in  Fig.  II- 3,  would  arrive  at  the  receiver  if  the  random  medium  in  the  right  half¬ 
space  were  removed  and  replaced  with  a  homogeneous  medium  similar  to  that  in  the  left  half¬ 
space.  Thus,  the  phase  fluctuation  could  be  obtained  if  the  equiphase  (or  isochron)  contours  for 
both  the  diffracted  wave  and  the  homogeneous  wave  were  known  for  each  event.  However,  the 
slowness  anomaly  data  for  each  event  can  be  used  only  to  obtain  the  gradient  of  the  equiphase 
contours  at  a  particular  point.  Thus,  it  is  necessary  to  obtain  the  equiphase  contours  from  the 
value  of  its  gradient  evaluated  at  nine  different  points  in  space. 

The  equiphase  contours  were  obtained  by  performing  a  least-squares  fit  of  their  gradients 
to  the  measured  gradients.  The  phase  fluctuation  for  the  j-th  event  and  at  the  i-th  subarray, 
denoted  S^,  j  =  1,  ...»  31,  i  =  1,  .  .  . ,  9,  is  then  obtained  from  the  equiphase  contours  in  the  usual 
manner.  The  estimate  of  the  transverse  autocorrelation  coefficient  for  the  phase  fluctuations, 
as  well  as  the  crosscorrelation  coefficient  between  the  amplitude  and  phase  fluctuations,  is  com¬ 
puted  in  a  manner  similar  to  that  for  the  amplitude  fluctuations.  In  addition,  the  estimate  for 
— 2  — 2 

B  and  S  ,  denoted  by  and  Pg,  respectively,  are  computed  in  the  standard  manner,  cf., 
Cramer.1  4 
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It  is  possible  to  plot  the  theoretical  relationship  between  B  /S  and  for  various  values 
of  the  wave  parameter,  as  shown  in  Fig.  II-5.  In  addition,  the  90 -percent-confidence  region  is 
shown  in  this  figure.  It  was  found  that  Pb  =  0.14,  Pg  =  0.27,  so  that  Pb/Pg  =  0.52,  and  that  the 
estimate  for  the  crosscorrelation  coefficient  is  0.2  3.  This  pair  of  measurements  is  shown  in 
Fig.  II-5,  along  with  the  90 -percent-confidence  region  for  D  =  6.  Thus,  the  wave  parameter  D 
was  estimated  to  be  equal  to  six,  as  this  appears  to  give  the  best  correspondence  between  the 
data  and  the  theory.  The  fact  that  the  measurement  lies  so  far  from  the  theoretical  curve  in 
Fig.  II -5  could  be  explained  solely  on  the  basis  of  the  finite  sample  sizes  used  in  the  estimation 
procedures.  However,  it  is  quite  likely  that  part  of  the  reason  for  the  departure  of  the  measure¬ 
ment  from  the  theoretical  results  is  due  to  some  of  the  basic  assumptions  not  being  satisfied. 

For  example,  the  assumptions  concerning  the  spatially  homogeneous  and  statistically  isotropic 
random  medium  are  clearly  not  satisfied.  In  view  of  this,  as  well  as  some  of  the  other  assump¬ 
tions  which  have  been  made,  the  agreement  of  the  measurements  with  the  theory  can  be  consid¬ 
ered  to  be  quite  good. 

The  results  obtained  for  the  estimates  of  the  transverse  autocorrelation  coefficients  of  am¬ 
plitude  and  phase  fluctuations  are  shown  in  Figs.  H-6(a)  and  -7(a),  respectively.  It  is  seen  from 
these  figures  that  the  coefficients  display  a  large  amount  of  variation,  so  that  some  smoothing 
was  employed  to  improve  the  stability  of  the  results,  as  shown  in  Figs.  II-6 (b)  and  -7(b).  Also 
depicted  in  these  figures  are  the  theoretical  results  given  by  Chernov^  for  Rb(f )  and  pre¬ 

sented  previously,  for  a  wave  parameter  D  equal  to  six  and  correlation  distance  a  of  1  2  km. 
There  seems  to  be  reasonably  good  agreement  between  the  experimental  and  theoretical  results 
in  these  figures,  so  that  it  may  be  concluded  that  a  good  estimate  for  the  correlation  distance 
at  LASA  is  12  km. 

It  is  now  possible  to  compute  an  estimate  for  the  extent  of  the  random  medium  L.  The  mean 
velocity  of  the  random  medium  is  c  =8  km/sec,  and  the  frequency  f  =  0.8  Hz,  so  that  to  =  27r(0.8)  = 
5.02  rad/sec.  Hence,  the  wavenumber  of  the  homogeneous  wave  is,  k  =  a)/cQ  =  27r(0.8)/8  = 

0.63  rad/km,  and  the  wavelength  \  =  27r/k  =  cQ/f  =  10  km.  Hence,  using  the  defining  equation 
for  D  we  have 

L  =  ka2D/4 

=  (0.63)(12)2(6)/4 
=  1  36  km 

It  is  also  possible  to  compute  an  estimate  for  the  standard  deviation  of  the  index  of  refraction, 

~~T  \  h  2  2 

(p  )  '  ,  using  the  two  relationships  given  by  Chernov  (pp.  82-83)  for  B  ,  S  .  Each  equation 

yields  a  different  value,  so  that  the  procedure  used  was  to  average  the  two  values  to  obtain  the 

2l/2. 

final  estimate.  Thus,  it  is  obtained  that  an  estimate  of  ( p  )  '  is  0.019.  Hence,  there  appear 
to  be  variations  in  the  index  of  refraction,  on  the  order  of  1.9  percent,  distributed  throughout 
the  crust  and  upper  mantle,  to  a  depth  of  more  than  100  km,  as  shown  in  Fig.  II—  8,  which  appear 
to  be  causing  the  P-wave  slowness  anomalies  between  subarrays  at  LASA.  What  can  be  inferred 
from  these  results?  It  must  be  emphasized  that  they  in  no  sense  prove  the  existence  of  a  random 
medium  under  LASA.  It  may  be  possible  to  find  a  set  of  corrugated  interfaces  which  can  explain 
the  amplitude  and  phase  anomalies.  This  has  not  been  possible  for  a  single  interface  (see 
Sec.I-F).  What  has  been  presented  is  an  estimate  of  random  media  parameters  which  are  ap¬ 
propriate  if  we  choose  to  describe  the  data  in  terms  of  Chernov  theory.  It  should  be  mentioned 
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that  similar  results  have  been  obtained  recently  by  Aki.  A  more  detailed  description  of  the 
present  results  will  be  presented  in  a  forthcoming  publication.1  ^ 

J.  Capon 

C.  LIMITATIONS  IMPOSED  BY  A  RANDOM  MEDIUM  ON  THE  ABILITY 
TO  MEASURE  dT/dA  WITH  LASA 

The  LASA,  as  well  as  other  array,  data  have  been  used  in  the  past  in  attempts  to  determine 
various  regions  of  heterogeneity  within  the  earth.  A  discussion  will  now  be  presented  of  the 
capability  of  the  LASA  to  determine  such  regions.  In  particular,  we  consider  only  those  methods 
based  on  the  analysis  of  the  P-wave  travel-time  derivative,  dT/dA.  In  these  investigations  dT/dA 
vs  A  is  measured  for  a  particular  set  of  azimuths  and  compared  with  a  standard,  such  as  the 
Jeffreys-Bullen  (J-B)  curve,  which  is  considered  reasonably  valid  for  all  azimuths.  Any  signi¬ 
ficant  differences  between  the  two  curves  is  interpreted  as  implying  the  existence  of  inhomoge- 
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neities  in  the  earth,  usually  deeper  in  the  mantle,  c.f.,  Chinnery  and  Toksoz,  Johnson, 
Corbishley.1  ^ 

Consider  the  various  sources  of  inhomogeneity.  There  may  be  various  inhomogeneous  re¬ 
gions  within  the  earth  that  cause  the  curve  for  dT/dA  vs  A  to  be  azimuthally  dependent,  when 
measured  at  a  particular  array.  In  addition,  this  curve  will  be  different  when  measured  at  dif¬ 
ferent  arrays.  These  inhomogeneities  may  be  described  broadly  as  due  to  the  following  causes. 

(1)  Random  part  of  local  crust  and  upper  mantle  structure  at  the  receiver. 

(2)  Deterministic  part  of  local  crustal  structure  at  the  receiver. 

(3)  Propagation  path  effects,  primarily  deeper  in  the  lower  mantle. 

(4)  Source  region  effects,  due  to  structures  such  as  dipping  plates  of  lithosphere. 

The  description  of  the  crust  and  upper  mantle  structure  at  LASA  as  a  random  medium,  given 

in  item  (1),  has  been  presented  in  Sec.  II-B,  'along  with  a  detailed  description  of  its  properties. 

This  represents  the  first  time  that  such  a  characterization  has  been  attempted  at  LASA.  The 

effects  due  to  the  structure  in  item  (2)  are  approximately  the  same  as  those  due  to  a  plane-dipping 

interface.  The  inhomogeneities  in  the  lower  mantle,  given  in  item  (3),  have  been  discussed  by 
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Chinnery  and  Toksoz  and  Johnson,  while  those  in  the  source  region,  given  in  item  (4),  were 

discussed  by  Davies  and  McKenzie20  and  Jacob.21  The  only  structure  listed  in  items  (1)  -  (4) 

which  can  cause  amplitude  and  phase  variations  between  subarrays,  for  a  given  event,  is  the 

random  medium  given  in  item  (1).  This  is  readily  apparent  and  may  be  seen  quite  clearly  from 

22 

some  of  the  ray  tracing  results  in  arbitrarily  heterogeneous  media  given  by  Julian.  These 

results  apply  quite  readily  to  the  structures  listed  in  items  (2),  (3)  and  (4),  and  show  that  these 

structures  must  produce  amplitude  and  phase  variations  that  are  constant  over  the  aperture  of 

\ 

LASA,  for  a  given  event.  It  has  been  shown  in  the  last  SATS  (Fig.  Ill- 3)  that  the  standard  devi¬ 
ations  of  these  variations  are  comparable,  or  larger,  than  their  average  for  each  event,  which 
must,  of  course,  be  due  to  the  structures  in  items  (2),  (3)  and  (4).  Thus,  it  follows  that  the 
structure  in  item  (1)  could  be  one  of  the  primary  causes  of  the  P-wave  anomalies  at  LASA.  How¬ 
ever,  contributions  to  the  P-wave  anomalies,  as  observed  on  individual  sensors,  due  to  the 
structures  listed  in  (2),  (3)  and  (4),  are  also  significant.  In  other  words,  all  of  these  structures 
may  be  of  importance  in  contributing  to  P-wave  anomalies  observed  on  individual  seismometers. 

A  method  which  has  been  used  to  measure  dT/dA  for  the  entire  LASA  is  to  note  the  P-wave 
arrival  times  at  the  center  of  the  subarrays  located  in  the  three  outer  rings,  D,  E  and  F,  and 
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to  find  the  slowness  of  a  plane  wave  whose  arrival  times  best  fit  the  measurements  in  a  least- 

squares  sense.  We  now  consider  the  effect  of  the  errors  introduced  by  a  random  medium  alone 

on  this  measurement.  The  standard  deviation  of  the  relative  error  in  dT/dA  for  this  method  has 
23 

been  given  by  Kelly  as 

/IV 

as  ■  V  N  t 

where  N  is  the  number  of  subarrays  used,  is  the  standard  deviation  of  the  arrival  time  meas¬ 
urement  (in  sec),  v  is  the  horizontal  phase  velocity  of  the  P-wave  for  a  particular  distance  (in 
km/sec),  and  i  is  the  aperture  of  the  array  (in  km).  It  can  be  shown  that  the  effect  of  the  random 
medium  is  to  introduce  timing  errors  whose  standard  deviation  is  about  0.10  sec.  Thus,  using 
v  =  20  km/sec,  corresponding  to  a  distance  of  about  80°,  N  =  12  and  H  =  200  km,  we  obtain 

_  nr  (o.i)  (20) 

s  ~  \l  12  200 

=  0.007 

Hence,  the  95-percent-confidence  limits,  due  to  the  effect  of  the  random  medium,  for  the  meas¬ 
urement  of  dT/dA  at  LASA  are  about  ±1.4  percent. 

Thus,  if  random  medium  results  are  applicable,  inferences  from  dT/dA  data  may  not  be 
significant  to  better  than  ±1.4  percent.  In  order  to  be  able  to  say  that  such  differences  are  due 
to  velocity  variations  in  the  lower  mantle,  item  (3),  or  in  the  source  region,  item  (4),  it  is  nec¬ 
essary  to  have  differences  that  exceed  ±1.4  percent,  and  the  effect  of  the  local  crustal  structure 

in  item  (2)  must  be  adequately  taken  into  account.  However,  this  level  of  the  difference  was  ex- 
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ceeded  in  the  analysis  of  Chinnery  and  Toksoz,  so  that  it  is  probably  that  the  differences  they 

observed  were  due  to  velocity  variations  in  the  lower  mantle,  as  they  claim,  and  not  due  to  the 

1  8 

random  medium.  It  is  also  possible  to  say  that  this  same  statement  applies  to  Johnson,  since 
he  used  data  at  the  Tonto  Forest  Observatory  (TFO)  and  the  effect  of  the  random  medium  may 
be  comparable  to  that  at  LASA. 

The  random  medium  has  a  small  effect  on  the  ability  to  measure  dT/dA  with  LASA,  and  only 
contributes  what  may  be  termed  a  small  experimental  error.  However,  this  may  not  be  true 
if  the  aperture  of  the  array  used  in  the  measurement  of  dT/dA  is  significantly  less  than  the  200-km 
aperture  of  LASA.  Some  care  must  be  exercised  in  ascribing  the  cause  of  differences  between 
the  measured,  and  the  J-B,  curve  for  dT/dA,  to  any  one  of  the  structures  listed  in  (2),  (3)  or  (4). 

J.  Capon 


D.  BODY-WAVE  MULTIPATHING 

The  two  most  fundamental  requirements  of  an  array  of  N  sensors  are  that  it  should  be  capable 
of  improving  the  SNR  by  something  approaching  \fN~ and  that  it  should  have  a  directional  capability 
to  discern  signals  arriving  in  different  directions  simultaneously.  Much  documentation  exists 
for  LASA's  ability  to  meet  these  requirements,  but  a  new  and  instructive  example  has  arisen  re¬ 
cently  which  is  particularly  relevant  to  the  problem  of  detecting  the  body  waves  of  events  in  the 
presence  of  interfering  (and  probably  larger)  signals. 

A  study  has  been  made  of  directions  of  approach  of  P-waves  at  individual  subarrays  of  LASA. 

1 

Earliest  results  from  it  were  reported  by  Capon  in  the  last  SATS.  More  recent  results  are  re¬ 
ported  elsewhere  in  Sec.II-A.  We  selected  a  few  events  for  further  study,  particularly  for 
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frequency  dependence  of  this  approach  angle.  In  most  cases  we  found  no  systematic  frequency 
dependence,  but  one  event  was  exceptional.  Observations  of  an  intermediate  focus  earthquake 
in  North  Colombia  (11  September  1966)  which  would  have  a  theoretical  dT/dA  of  7.6  sec/deg  and 
an  azimuth  of  134°,  showed  that  energy  at  the  low  frequency  end  of  the  spectrum  (0.6  to  0.8  Hz) 
consistently  approached  from  a  direction  significantly  different  from  that  of  the  high  frequency 
(1.1  to  1.3  Hz)  energy.  In  Fig.  II-9  we  show  a  high  resolution  plot  of  the  frequency  wave-number 
spectrum  at  0.8  and  1.2  Hz  for  subarray  B4.  Neither  frequency  approaches  at  the  theoretical 
direction,  but  the  absolute  position  of  individual  power  peaks  means  little,  as  geological  structure 
in  the  vicinity  of  a  subarray  may  severely  bias  the  arrival  direction.  However,  differential  posi¬ 
tions  should  be  significant.  Thus,  we  have  plotted  for  each  subarray  low  and  high  frequency 
points  in  an  ’’array  diagram.”  The  vectors  joining  mean  low  frequency  and  mean  high  frequency 
positions  for  all  9  subarrays  within  the  C-ring  are  shown  in  Fig.  II- 1 0,  where  it  is  apparent  that 
they  do  not  have  a  zero  mean.  The  effect  is  dominantly  an  azimuthal  one.  (We  infer  that  low  and 
high  frequency  energy  approach  LASA  on  azimuths  differing  by  10°  to  20°.  )  We  cannot  tell  whether 
either  direction  of  approach  is  the  ’’correct”  one  without  a  detailed  analysis  of  subarray  structure, 
but  for  the  purpose  of  this  study  this  is  unnecessary.  Both  signals  arrive  at  almost  the  same  time 
so  the  most  reasonable  explanation  would  be  that  both  follow  non-great  circle  paths,  one  each  side 
of  the  theoretical  path.  The  spectra  of  these  two  signals  are  significantly  different.  Propagation 
in  regions  of  strong  lateral  velocity  gradients  could  account  for  this,  but  this  is  speculation. 

We  have  also  analyzed  the  center  seismometers  of  subarrays  out  to  the  C-ring  and  these 
taken  as  a  whole  also  show  the  same  effect.  Figures  II- 11  (a)  and  (b)  show  the  conventional 
frequency-wave  number  plots  at  0.6  and  1.2  Hz. 

We  believe  that  the  source  of  this  multipathing  has  to  be  sought  beyond  the  immediate  sub¬ 
structure  of  the  array,  as  the  effect  does  not  appear  on  events  from  other  regions.  We  see  pos¬ 
sible  causes  as  being  restrained  to  be  in  the  middle  mantle  where  rays  are  traveling  nearly  hori¬ 
zontally.  Lateral  velocity  gradients  at  depths  of  up  to  1200  km  underneath  Mexico,  the  Gulf  of 
Mexico  and  Florida  might  account  for  two  refracted  paths.  Whatever  the  cause,  it  strongly 
affects  coherence  of  the  P-wave  across  the  array  and  leads  to  the  need  for  a  certain  care  in 
describing  array  directional  capability  when  the  primary  signal,  at  different  frequencies,  is 

not  from  one  unique  direction.  ~  _ 

^  D.  Davies 

J.  Capon 


E.  VELOCITY  VARIATIONS  IN  THE  LOWER  MANTLE 

1 

The  June  1972  SATS  contained  a  description  of  a  P-wave  travel  time  experiment  which  was 
conducted  using  data  from  deep  focus  earthquakes  exclusively.  One  of  the  most  significant  con¬ 
clusions  was  that  data  at  distances  beyond  70°  or  so  were  found  to  contain  anomalies  which  could 
be  explained  satisfactorily  only  in  terms  of  lateral  velocity  variations  in  the  deep  mantle.  These 
data  have  now  been  analyzed  in  more  detail,  and  a  refined  estimate  of  the  nature  of  the  variations 
has  been  obtained. 

Owing  to  the  uneven  distribution  of  seismological  observatories  and  deep  earthquakes,  the 
sampling  of  the  mantle  provided  by  these  data  is  also  quite  uneven,  and  it  is  impossible  to  deter¬ 
mine  uniquely  the  complete  three-dimensional  velocity  structure  of  the  mantle.  What  can  be 
determined  is  the  average  travel  time  residual  for  each  of  a  number  of  "bundles”  of  rays  following 
nearly  identical  paths  from  a  seismic  region  to  a  group  of  stations.  For  each  path,  a  Student's 
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t  test  was  used  to  evaluate  the  hypothesis  that  the  travel  time  (after  station  corrections  were 
applied)  is  the  value  given  by  the  overall  travel  time  curve  and  that  deviations  from  this  curve 
can  be  attributed  to  random  measuring  errors.  Those  ray  paths  for  which  the  hypothesis  could 
be  rejected  at  the  99.5  percent  confidence  level  were  found  almost  exclusively  at  distances  beyond 
70 implying  that  velocity  variations  near  the  bottoming  points  of  the  rays  produce  the  travel 
time  anomalies.  Other  effects,  such  as  near-source  or  near-receiver  velocity  anomalies  or 
event  mislocation,  are  incapable  of  explaining  the  pronounced  distance  dependence  observed. 

The  actual  details  of  the  velocity  distribution  cannot,  however,  be  determined  precisely, 
because  rays  travel  a  great  distance  at  approximately  the  same  depth  near  their  turning  points. 
Figure  11-12  shows  the  regions  of  the  lower  mantle  sampled  by  the  various  ray  bundles  (defined 
arbitrarily  by  the  middle  30°  of  the  rays)  and  indicates  which  paths  correspond  to  early  and  late 
arrivals.  Several  interesting  features  of  this  figure  should  be  noted.  First,  where  regions 
overlap  on  the  figure,  they  are  generally  consistent  with  each  other  (e.g.,  regions  1-6  and  1-7, 
and  regions  2-2  and  2-3).  This  consistency  is  encouraging,  in  that  it  supports  our  argument  that 
the  travel  time  anomalies  originate  in  the  deep  mantle  and  are  not  the  result  of  some  other  type 
of  systematic  error.  An  apparent  inconsistency  exists  between  regions  1-3  and  1-16,  but  this 
is  not  surprising  in  view  of  the  uncertainty  as  to  exactly  where  the  travel  time  anomalies  actually 
originate.  In  fact,  the  rays  following  path  1-3  also  pass  through  regions  2-2  and  2-3  further  to 
the  north,  and  it  is  likely  that  the  travel  time  anomaly  actually  originates  there.  Another  inter¬ 
esting  feature  of  Fig.  11-12  is  a  correlation  between  the  anomalies  in  the  two  distance  ranges 
(e.g.,  regions  1-7  and  2-7,  regions  1-4  and  2-4,  and  regions  1-5  and  2-6),  suggesting  that  the 
structure  in  the  deep  mantle  has  a  spatial  11  coherence1'  of  at  least  a  few  hundred  km,  vertically. 

The  mean  travel  times  for  the  various  ray  bundles  show  a  variation  of  about  1.5  sec  for  rays 
bottoming  below  2600  km.  This  number  is,  of  course,  subject  to  some  uncertainty,  but  it  is 
nearly  certain  that  the  true  travel  time  varies  by  at  least  1  sec.  The  amount  by  which  the  actual 
velocity  in  the  mantle  varies  depends  on  the  size  of  the  regions  within  which  the  variations  occur. 
Figure  11-12  suggests  that  the  size  of  some  of  the  anomalies,  at  least,  is  about  1000  km  or  less, 
in  which  case  the  velocity  must  vary  by  at  least  1  percent.  This  figure  is  a  lower  bound,  both 
because  we  have  probably  overestimated  the  scale  of  the  inhomogeneities  and  because  we  are 
measuring  averages  of  the  velocity  in  rather  large  regions  and  some  cancellation  of  the  effects 
of  positive  and  negative  velocity  anomalies  is  to  be  expected.  Combined  interpretation  of  travel 
time  and  dT/dA  measurements  can  probably  improve  the  resolution  of  structural  details. 

B.  R.  Julian 

M.  K.  Sengupta 

F.  EFFECTS  OF  THE  SEDIMENTS  AT  LASA  ON  SHORT-PERIOD  TELESEISMS 

Data  showing  that  the  boundary  between  the  sediments  and  basement  under  LASA  is  a  strong 
scatterer  of  S-waves  were  presented  in  the  June  1972  SATS*  It  was  further  indicated  that  the 
boundary  was  dipping  and  that  the  resulting  structure  might  account  for  some  of  the  short  period 
waveform  anomalies  observed  at  the  array.  Further  data  have  now  been  analyzed  in  order  to 
detail  the  nature  of  the  boundary  so  as  to  calculate  its  effect  on  emerging  P-waves. 

Figure  II-l  3  shows  the  azimuthal  dependence  of  the  amplitude  ratio  of  the  scattered  S-wave 
to  the  scattering  P-wave.  The  pattern  indicates  that  the  boundary  is  essentially  planar  and 
dipping  to  the  southeast.  To  determine  the  dip,  the  same  ratio  was  measured  for  events  in  the 
down-dip  direction  at  various  distances.  The  data  are  shown  in  Fig.  11-14.  The  theoretical 
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curves  for  the  structure  shown  in  the  figure  were  computed  using  the  Jeffreys-Bullen  dT/dA  vs  A 

values  and  plane  wave-plane  interface  reflection  coefficients,  an  example  of  which  is  shown  in 

Fig.  11-15.  The  S-wave  velocity  in  the  sediments  was  determined  by  using  the  P-scattered  S 
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arrival  times  and  the  P-wave  velocity  and  sediment  thickness  given  by  Brown  and  Poort  from 
well  log  data.  The  calculated  S-wave  velocity  indicates  a  Poisson’s  ratio  of  0.38.  Such  a  value 
is  indicative  of  the  clastic  sediments  at  LASA.  The  best  fitting  curve  in  Fig.  11-14  is  for  a  dip 
of  10  degrees.  The  technique  is  in  fact  quite  sensitive  to  the  dip  of  the  layer. 

To  examine  the  effects  of  the  structure,  the  amplitude  and  direction  of  the  primary  wave  and 
higher  order  multiply  reflected  waves  were  calculated.  Two  examples  are  shown  in  Figs.  II - 1 6 
and  -17.  For  events  whose  P-waves  emerge  from  the  down-dip  direction,  a  first  multiple  P-wave 
with  an  amplitude  of  about  0.4  of  the  primary  wave  (it  will  be  out  of  phase)  will  be  mixed  with  the 
primary  at  an  offset  equal  to  the  two  way  travel  time  in  the  crust.  Higher  order  multiples  will 
be  down  by  an  order  of  magnitude  as  will  the  first  order  multiple  for  events  whose  P-waves 
emerge  from  the  up-dip  direction.  The  two-way  travel  time  for  a  P-wave  in  the  crust  will  be 
approximately  2.2  5  sec.  Figure  11-18  shows  the  P  and  SV  components  of  2  events  from  the  down- 
dip  direction.  The  SV  components  have  been  shifted  so  that  they  are  directly  beneath  the  scatter¬ 
ing  P-wave.  P-wave  motion  without  corresponding  S-wave  motion  indicates  that  the  P-wave 
motion  was  generated  above  the  sediment-basement  boundary  and  is  not  source  motion.  The 
arrivals  marked  Pr  have  the  corresponding  time  amplitude  and  phase  of  the  proposed  first  order 

multiple.  Allowing  for  the  fact  that  the  depth  of  the  boundary  varies  from  1800  to  3400  m  over 
24 

the  array  and  that  its  dip  and  direction  is  probably  nearly  random  from  subarray  to  subarray 

24 

with  magnitudes  of  at  least  10  degrees,  the  observed  short-period  waveform  anomalies  may 
be  due  to  the  influence  of  the  layer.  However,  only  about  30  percent  of  the  amplitude  anomalies 
at  LASA  can  be  explained  by  the  structure.  Thus,  while  the  upper-crust  has  a  significant  effect 
on  short-period  waves,  it  is  only  partly  responsible  for  the  behavior  of  short-period  P-waves 
observed  at  the  array.  x.  Landers 
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Fig.  II— 1 .  Subarray  diagrams  for  B-ring  at  LASA:  (a)  site  Bl;  (b)  site  B2; 
(c)  site  B3;  (d)  site  B4.  Outer  circle  radius  is  10  sec/deg. 
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Fig.  II-2.  Subarray  average  (a),  spatial  average  (b),  and  theoretical  array  diagram  (c) 
for  a  plane-dipping  interface.  Outer  circle  radius  is  10  sec/deg.  For  theoretical  array 
diagram  (c),  dip  =  11  °,  dip  direction  =  316°. 


MEDIUM  MEDIUM 

(lower  mantle)  (crust  and  upper  mantle) 


Fig.  II—  3.  Illustration  of  model  used  by  Chernov  for  wave  propagation 
in  a  random  medium. 
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Fig.  II-4.  Actual  and  simplified  compressional  velocity  structures  in  the  crust 
and  upper  mantle. 
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Fig.  II- 5.  Theoretical  relationship  between  amplitude -relative-to-phase 
fluctuation  power  ratio,  and  crosscorrelation  coefficient  between  ampli¬ 
tude  and  phase  fluctuation,  and  the  90  percent  confidence  limits  of  the 
measurement. 
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Fig.  II-6.  Comparison  of  unsmoothed  and  smoothed  estimates  with  theoretical 
curve  for  transverse  autocorrelation  coefficient  of  amplitude  fluctuations.  Wave 
parameter  =  6;  correlation  distance  =  12  km. 
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Fig.  II- 7.  Comparison  of  unsmoothed  and  smoothed  estimates  with  theoretical 
curve  for  transverse  autocorrelation  coefficient  of  phase  fluctuations.  Wave 
parameter  =  6;  correlation  distance  =  12  km. 


49 


COM  PRESS  I  ON  AL  VELOCITY 
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Fig.  II- 8.  Diagrammatic  representation  of  crust  and  upper  mantle 
structure  under  LASA. 
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Fig.  II-9.  High  resolution  frequency-wavenumber  spectra  at  0.8  (a)  and  1.2  Hz  (b) 
for  subarray  B4  for  a  North  Colombian  event.  Power  peak  is  at  61.8,  and  75.2  dB, 
respectively,  above  one  quantum  unit  level.  Contours  are  2  dB  apart. 
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Fig.  11-10.  Distribution  of  vectors  joining  low  to  high  frequency  directions 
of  approach  in  an  array  diagram  for  all  subarrays  within  the  C-ring.  Mean 
implies  an  azimuthal  separation  of  about  20°. 
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(a) 


(b) 

Fig.  H-ll.  Conventional  frequency-wave  number  spectra  at  0.6  (a)  and  1.2  Hz  (b) 
for  center  seismometers  of  all  subarrays  out  to  the  C-ring.  Peaks  are  at  (10.2, 
131°)  and  (9.1,142°),  respectively.  Intersection  of  circle  and  line  is  expected 
wavenumber. 
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LATE  EARLY 


LATE  EARLY 


Fig.  11-12.  Regions  of  deep  mantle  sampled  by  "bundles"  of  observed 
rays,  (a)  85°  <  A  <  100°  (depth  2600-2900  km);  (b)  70°  <  A  <  85°  (depth 
2000-2600  km). 
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Fig.  II -1  3.  Azimuthal  dependence  of  scattered 
S-wave  amplitude  indicating  a  layer  dipping  to 
the  southeast. 


AZIMUTH  *  140°  TO  160° 


Fig.  11-14.  Distance  dependence  of  scattered  S-wave  amplitude 
indicating  a  dip  of  10  degrees  exists  at  the  scattering  boundary. 
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Fig.  11-15.  Amplitude  of  transmitted  P-  and  S-waves  and  reflected  P- 
and  S-waves  as  a  function  of  incidence  angle  for  the  structure  shown. 


Fig.  11-16.  Ray  paths  and  amplitudes  of  scattered  waves  for  events 
emerging  from  down-dip  direction. 
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Fig.  11-17.  Same  as  Fig.  11-16  for  events  emerging  from  up-dip  direction. 
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A  =  73°  AZIMUTH  =  145° 


T  +  1.3  sec 


Fig.  11-18.  Seismograms  of  P- waves  from  events  in  down- dip  direction 
showing  first  order  multiple  (P  ). 
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III.  SOURCE  STUDIES 


A.  FAULT  PLANE  SOLUTIONS  IN  ASIA 

Because  of  the  failure  of  certain  Asian  events  (assumed  to  be  earthquakes)  to  discriminate 
from  explosions  through  the  mb:Mg  ratio  (T.  Landers  in  Ref.  l),we  have  begun  a  study  into  the 
seismicity  and  tectonics  of  central  and  eastern  Asia.  In  the  initial  phase  of  this  study,  we  have 
collected  first  motion  data  in  an  effort  to  determine  the  orientation  of  various  tectonic  zones 
through  fault  plane  solutions.  Since  no  clear,  unique  pattern  has  arisen  to  date,  we  shall  only 
review  our  data  and  suggest  tentative  trends. 

The  seismicity  map  of  the  NOS,  the  base  map  of  Figs.  Ill- 1  and  -2,  shows  the  location  of 
earthquakes  for  the  1962-1 969  period.  From  this  map  it  is  clear  that  the  seismicity  of  China, 
Tibet,  Sinkiang,  and  the  Mongolias  shows  no  well  defined  trends.  Linear  zones  which  charac¬ 
terize  plate  boundaries  in  other  regions  are  not  immediately  evident  here.  Rather  the  epicenters 
are  widely  scattered,  except  in  a  few  regions  where  dense  pockets  of  activity  occur  (37 °N,  115°E 
near  Peking,  37°N,  95°E  in  Tsinghai  province,  and  30°N,  95°E  in  eastern  Tibet).  The  historical 
record^  shows  a  linear  concentration  (the  Lanchow-Kunming  Fault  zone)  of  seismicity  along  the 
103°E  meridian,  beginning  at  24°N  and  running  to  34°N  where  it  disperses.  One  branch  of  this 
activity  then  runs  northeastward  in  an  S-shaped  pattern  ending  near  Peking.  The  fact  that  such 
patterns  show  only  when  historical  data  are  considered  indicates  that  the  seismicity  of  these  re¬ 
gions  is  a  relatively  low  level  when  compared  with  the  Hindu-Kush.  Further,  the  lack  of  earth¬ 
quakes  large  enough  for  reliable  fault-plane  solutions  at  teleseismic  distances  makes  the  prob¬ 
lem  of  tectonic  interpretation  more  acute. 

Using  the  first  motions  recorded  on  long  period  instruments  of  the  WWSSN,  we  have  con¬ 
structed  fault  plane  solutions  of  some  18  earthquakes.  The  solutions,  along  with  nine  reported 
3  4 

by  Fitch  and  five  by  Nowroozi,  are  summarized  in  Fig.  Ill— 1 .  The  circles  represent  the  lower 
hemispheres  of  equal  area  projections  of  the  first  motion  data,  the  solid  areas  represent  areas 
of  compression.  The  heavy  lines  on  this  map  represent  "main  deep  faults”5  either  traced  (solid) 
or  hypothetical  (dashed). 

In  some  cases  these  fault  lines  are  parallel  to  possible  fault  planes  suggested  by  the  first 
motion  data.  Earthquakes  numbered  1,  2,  9,  14,  22,  23  and  F6  are  such  examples.  Other  solu¬ 
tions,  however,  such  as  17  and  FI,  show  fault  planes  nearly  normal  to  the  mapped  faults.  Still 
others,  notably  15  and  16,  do  not  lie  near  any  prominent  surface  fault  feature.  As  with  the  seis¬ 
micity,  these  fault  plane  solutions  yield  no  clear,  simple  pattern  for  the  tectonics  of  this  Asian 
region.  Most  of  Fitch's  solution  along  the  Himalayan  front  and  Nowroozi' s  in  the  Hindu  Kush 
indicate  thrust  faulting  with  compression  from  the  south  and  west.  The  solutions  falling  in  north¬ 
ern  Burma,  eastern  Tibet  and  southwestern  China  are  quite  a  mixed  bag.  Normal,  thrust  and 
strike-slip  mechanisms  are  all  present  and  we  cannot,  for  discrimination  purposes,  exclude  the 
possibility  of  any  of  these  types.  Although  Wilson^  has  suggested  that  the  Lanchow-Kunming 
fault  zone  represents  a  plate  boundary,  our  solutions  in  or  near  this  zone  do  not  indicate  what 
the  relative  motions  should  be. 

We  have  attempted  to  construe  a  motion  pattern  from  poles  of  the  auxiliary  planes  of  the  so¬ 
lution  in  Fig.  Ill— 1 .  We  start  with  a  solution  that  shows  a  fault  plane  striking  parallel  to  a  mapped 
fault  and  then,  in  cases  where  neither  plane  is  the  obvious  fault  plane,  pick  that  plane  whose 
motion  vector  is  most  compatible  with  surrounding  motion  vectors.  Of  course,  such  an  exercise 
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does  not  yield  a  unique  set  of  motions  (or  relative  slip  vectors).  The  motion  vectors,  which  are 
pointing  northerly  near  the  bottom  center  of  the  map,  swing  to  the  west  in  the  center  of  the  map. 
These  vectors  are  roughly  parallel  to  the  structural  geologic  trends  of  the  region  around  the 
high  angle  intersection  of  the  arcs  representing  the  Himalayan  and  Burmese  fronts.  It  appears 
as  though  this  portion  of  Asia  was  being  bent  (or  stretched)  about  a  point  in  Assam. 

Otherwise,  the  interpretation  of  Fig.  Ill— Z  indicates  that  northern  China  is  subject  to  east- 
west,  rather  than  north-south,  shear.  In  these  isolated  fault  plane  solutions  this  is,  at  present, 
an  arbitrary  interpretation.  More  specifically,  we  feel  we  can  state  that  the  region  of  Tibet 
which  yields  explosion-like  earthquakes  is  one  of  complicated  tectonics  not  dominated  by  one 
type  of  focal  mechanism.  j  pQson 

S.  Das 


B.  A  THRUST  EVENT  IN  SOUTHERN  IRAN 

On  10  April  1972,  a  large  destructive  earthquake  (NOS  Mg  6.9)  occurred  near  the  town  of 
Ghir,  some  150  km  south  of  Shiraz,  Iran.  The  main  shock  was  followed  by  many  aftershocks, 
of  which  25  were  reported  during  the  period  10  April  to  25  May  in  the  NOS  Earthquake  Data 
Reports. 

Figure  III— 3  shows  the  fault -plane  solution  determined  from  the  polarities  of  initial  P-wave 
motions  recorded  atWWSSN  stations.  It  proved  impossible  to  draw  two  mutually  orthogonal  planes 
through  the  data  when  the  focus  was  set  in  the  upper  mantle  (P-wave  velocity  =8.0  km-sec  ), 

_  j 

thus  it  was  placed  in  the  crust  (6.8  km-sec  ).  This  has  the  effect  of  plotting  all  stations  closer 
to  the  center  of  the  equal  area  projection,  because  of  refraction  at  the  Moho.  It  was  then  pos¬ 
sible  to  draw  two  nodal  planes  through  the  data.  Only  two  dilatations  were  recorded  and  a  thrust- 
type  event  is  clearly  indicated.  The  upper  plane  (I)  is  well  constrained:  the  lower  one  (II)  is 
not.  It  can  take  any  position  between  IIA  and  IIB  and  still  satisfy  both  the  data  and  the  orthog¬ 
onality  condition. 

The  seismic  moment  was  calculated  from  the  ratio  of  observed  to  theoretical  spectral  am¬ 
plitudes  of  the  long-period  vertical  Rayleigh  waves  and  the  mean  value  obtained  was  (1.5  ±  0.4)  X 
102^  dyne-cm.  Theoretical  spectra  were  calculated  using  the  method  of  Saito  in  a  flat  continental 
earth  model,  with  source  depth  assumed  to  be  30  km  and  source  time  dependence  assumed  to  be 
that  of  a  step  function. 

The  main  shock  was  relocated  (28.51  °N,  52.82°E)  and  the  residuals  used  as  station  correc¬ 
tions  in  relocating  the  aftershocks.  This  proved  to  be  remarkably  successful  in  restricting  these 

to  a  fairly  narrow  band  (Fig.  III-4)  in  comparison  with  the  NOS  locations,  which  were  spread  over 
2 

an  area  of  5000  km  .  On  the  same  diagram  are  shown  the  observed  fault  breaks  on  the  surface. 

The  trend  of  these  agrees  well  with  those  of  the  aftershock  distribution.  There  is  no  indication 

as  to  which  of  the  nodal  planes  is  the  fault  plane.  To  agree  with  the  fault  break/aftershock  trend 

it  could  be  either  plane  I  or  plane  IIB,  in  which  case  the  event  has  no  strike-slip  component. 

The  length  of  the  aftershock  distribution  (85  km)  is  taken  as  that  of  the  fault  plane  and  the 

depth  assumed  (no  depth  phases  were  reported)  to  be  30  km.  This  gives  the  area  of  the  fault 

plane  A  as  (85  X  30)/cos  6  km  ,  where  6  is  the  dip  of  the  (unresolved)  fault  plane  which  can  be 

11  -2 

either  54°  (IIB)  or  36°  (I).  Taking  |jl  =  5x10  dyne-cm  ,  we  have,  since  the  seismic  moment 
26 

Mq  =  puA  =  1.5x10  dyne-cm  (where  \i  is  the  shear  modulus  at  the  source  depth),  that  the 
minimum  value  of  mean  relative  displacement  u  along  the  fault  plane  is  (12  cos  54°)  =  7  cm. 
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Taking  Eg  =  5.4  +  2.8  as  a  measure  of  the  seismic  energy  release,  we  have,  from  m^  = 
6.1  and  E,  the  total  energy  release  =  cfuA.  Then  the  apparent  average  stress,  the  product  of 
mean  stress  on  the  fault  plane  and  the  seismic  efficiency  factor  is  10  bars. 

R.  North* 


C.  OBSERVATIONS  OF  DEPTH  PHASES 

The  phases  pP  and  sP  are  of  the  utmost  importance  in  event  location  and  identification. 

Depth  determination  by  times  of  P  arrivals  is  never  entirely  satisfactory.  Ln  particular,  for 

small  events,  where  less  than  10  stations  may  be  reporting,  the  depth  estimate  is  necessarily 

unreliable,  especially  for  shallow  events.  Under  these  circumstances  the  presence  of  pP  or  sP 

1 

is  most  useful.  In  the  last  SATS,  for  instance,  Landers  showed  the  impact  that  pP  made  on 
identifying  an  event  which  could  not  be  ascribed  a  depth  to  within  50  kms  by  P-waves  alone. 

Unfortunately  pP  and  sP  are  somewhat  wayward  phases.  I  do  not  have  extensive  statistics 
on  depth  phases,  but  based  on  limited  experience  it  appears  that  frequently  these  phases  are  not 
recorded  adequately  enough  to  permit  their  confident  identification.  This  may  be  because  the 
phase  occurs  in  a  high  coda  level,  but  I  present  here  some  evidence  to  show  that,  at  least  for 
one  event,  pP  and  sP  are  absent  in  the  short-period  band  in  a  case  where  the  coda  level  is  low 
and  the  focal  mechanism  gives  no  indication  of  any  unusual  circumstances.  The  event  was  a 
Hindu  Kush  earthquake  (m^  about  6.0)  of  20  January  1972.  The  locations  from  LASA  and  NORSAR 
were  36 °N,  7  0°E  and  34°N,  7  3°E,  respectively.  Neither  of  the  arrays  reported  a  depth  phase 
(Fig.  Ill — 5)  and  the  short-period  P-waves  were  fairly  simple.  It  could  be  argued  that  the  LASA 
beam  shows  a  signal  about  53  sec  after  P  and  that  this  is  in  fact  pP,  but  by  the  same  token  a 
signal  about  20  sec  after  P  could  be  called  pP.  In  addition,  surface  waves  were  very  poorly 
excited.  Thus,  on  the  face  of  it  this  event  needed  further  investigation,  especially  as  its  array 
locations  fell  relatively  close  to  the  borders  of  the  Soviet  Union  (Fig.  III-6). 

Hindsight,  in  the  form  of  the  NOS  bulletin,  tells  us  that  the  event  was  located  at  36.4°N, 

70.7 °E  and  was  213  km  deep,  which  explains  the  lack  of  surface  waves.  However,  before  a  net¬ 
work  depth  determination  became  available  there  was  no  short-period  body  wave  evidence  that 
the  event  was  not  at  the  surface.  In  Fig.  III-7  is  a  comparison  between  short-period  P  at  College, 
Alaska,  for  this  event  and  an  event  at  36.3°N,  70.4°E,  h  =  22  5  km  on  29  January  1968.  The  clear 
pP  for  the  latter  event  demonstrates  the  marked  variability  in  reflectivity  in  this  region.  The 
problem  in  this  case  is  not  insoluble.  Although  no  station  among  40  or  50  WWSSN  and  array  lo¬ 
cations  examined  gave  anything  approaching  a  credible  short-period  pP,  analysis  of  long-period 
records  has  revealed  evidence  for  long-period  sP,  but  nothing  very  conclusive  for  long-period 
pP.  In  Fig.  Ill— 8  are  6  long-period  records  showing  P  and  sP  and  the  limited  evidence  for  pP 
(at  Chiengmai,  Thailand).  Figure  III-9  is  a  LASA  beam  showing  rather  clearly  both  sP  and  sPP. 

It  appears  that  this  event  displays  an  unusual  combination  of  circumstances.  The  focal 
mechanism  is  an  overthrust  on  a  fault  plane  very  nearly  vertical  or  horizontal.  Long-period 
pP  to  stations  in  the  northern  hemisphere  is  rather  poorly  excited  being  near  a  nodal  plane. 
Short-period  pP  and  sP  are  rendered  unrecognizable,  either  because  of  strong  scattering  near 
the  Earth’s  surface  (topography  is  very  marked  in  this  region),  or  by  some  defocusing  owing  to 
lateral  variations  of  elastic  properties  near  source,  or  by  both.  Probably  both.  The  event  high¬ 
lights  the  necessity  to  take  a  very  broad  view  of  the  source. 

D.  Davies 

*  University  of  Cambridge. 
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D.  A  STOCHASTIC  EARTHQUAKE  SEQUENCE 

Knopoff7  has  set  down  a  stochastic  model  for  earthquake  occurrence  in  which  the  potential 
energy  of  deformation  is  taken  as  the  appropriate  state  variable.  The  concept  has  recently  been 

g 

tested,  with  favorable  results,  against  a  laboratory  device  built  to  model  a  seismic  region.  The 
stochastic  model  is  described  by  the  integrodifferential  equation 


where 


and 


*2P+ff  H 


=  fE] 

^E 


JMAX 


T(X,  E)  X2(X)  P(X)  dX 


(III-l) 


P(E)  dE  is  the  probability  that  the  system  is  in  a  state  of  deformation 
with  energy  between  E  and  E  +  dE, 

A^E)  dt  is  the  probability  that  an  event  occurred  in  the  time  interval 
dt  if  the  system  is  in  energy  state  E, 

T(X,  E)  dE  is  the  conditional  probability  that,  if  the  system  is  in  the 

pre-shock  state  of  deformation  X  and  a  shock  does  occur 
then  the  final  state  of  deformation  lies  between  E  and 
E  +  dE. 

a  is  the  rate  at  which  deformation  energy  is  added  to  the  sys¬ 
tem  and  is  taken  to  be  constant. 


I^mAX  is  the  maximum  energy  of  deformation  which  the  system  can 
sustain. 


Equation  (III— 1 )  is  based  on  the  assumption  that  the  only  way  an  energy  state  E  can  be  reached 
is  from  below  through  the  build-up  of  deformational  energy  or  from  above  through  the  occurrence 
of  an  earthquake  at  some  greater  energy  state.  In  Eq.  (Ill— 1 )  all  the  probabilities  are  taken  to  be 
stationary. 

The  verification  of  Eq.  (Ill— 1 )  in  nature  would  be  based  on  an  earthquake  catalog  of  an  auton¬ 
omous  seismic  region  where  the  rate  of  deformational  energy  increase  is  known.  The  latter  two 
conditions  may  not,  in  general,  be  frequently  met.  Tom  Simkin  (Smithsonian  Institution)  and  I 
have  recently  completed  a  rather  detailed  study  of  an  earthquake  swarm  which  occurred  during 
June  of  1968  in  the  Galapagos  Islands.  We  have  attributed  the  source  of  the  swarm  to  the  col¬ 
lapse  of  the  caldera  of  the  island  volcano  Fernandina.  From  our  seismicity  lists  we  have  com¬ 
puted  the  seismic  energy  release  as  a  function  of  time  using  surface  wave  magnitudes  (M  )  and 
9  s 

the  formula 

logE  =  11.8  +  1.5  M  .  ( III— 2 ) 

Figure  III— 1 0  summarizes  the  seismicity  of  the  Galapagos  swarm.  The  histogram  at  the  top  of 
the  figure  indicates  the  number  of  events  that  occurred  during  each  day,  while  the  cumulative 
energy  based  on  Eq.  (III-2)  is  plotted  in  the  lower  part  of  the  figure  on  both  logarithmic  and  linear 
scales.  It  is  clear  from  the  linear  energy  vs  time  diagram  that  the  rate  of  energy  release  is 
nearly  constant  over  about  a  seven  day  period.  We  suggest  that  gravitational  potential  energy  is 
supplied  to  the  caldera  block  by  the  withdrawal  of  a  supporting  magma  from  below  at  a  constant 
rate  and  that  this  phenomenon  is  reflected  in  the  constant  rate  of  seismic  energy  release.  If 
these  assertions  are  correct,  then  the  physical  conditions  of  autonomy  and  constant  rate  of  de¬ 
formational  energy  increase  necessary  for  the  verification  of  Eq.  (Ill— 1 )  are  met. 
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Following  Knopoff,  we  have  tested  Eq.  (III-l)  by  evaluating  the  left  and  right  hand  sides  at 

twenty  equal  energy  increments  from  E  =  0  to  E  =  over  the  period  13-19  June  1968.  We 

used  a  time  increment  of  1  min.  and,  from  Fig.  HI — 1 1  have  set  ot  =  0.6  x  10*^  ergs/min. 

Throughout  this  test  we  take  the  seismic  energy,  E  ,  to  represent  the  energy  of  deformation, 

s  20 

E,  through  a  seismic  efficiency  constant.  We  have  set  =  0*4  X  10  ergs,  the  seismic 

energy  of  the  largest  earthquake  of  the  sequence.  At  each  minute  of  the  seven  day  period  we 
have  computed  the  deformational  energy  of  the  system,  E(t),  using 

rt-l 

E(t)  =  (a  +  cvt)  -  \  Egdt  (III-3) 

20 

where  was  set  equal  to  0.15  X  10  ergs.  Using  our  list  of  origin  times  and  surface  wave 
magnitudes  for  341  earthquakes  during  the  period  and  Eqs.  (III-l  and  -2),  the  evaluation  of  both 
sides  of  Eq.  (III-l)  is  straightforward. 

In  Fig.  III-l  1  the  values  of  the  right-  and  left-hand  sides  of  Eq.  (III-l)  are  plotted  for  the 
20  energy  increments  used.  The  agreement  between  the  two  sets  of  values  is  good  and  Eq.  (III-l) 
is  taken  as  representative  of  the  statistical  aspects  of  the  swarm.  We  have  also  evaluated 
Eq.  (III-l)  over  sub-intervals  of  less  than  seven  days  and  found  that,  although  Eq.  (III-l)  is  still 
satisfied,  the  shape  of  the  probability  density  functions  P(E)  and  X^E)  clianges-  This  implies 
that  although  the  statistics  of  the  swarm  are  not  time  stationary  they  are,  over  any  time  interval, 
related  through  Eq.  (III-l).  That  the  swarm  is  nonstationary  can  easily  be  seen  in  Fig.  III-l 2 
where  we  show  the  detail  of  the  first  five  days  of  the  major  portion  of  the  swarm.  It  is  clear 
from  this  figure  that  the  swarm  is  not  stationary,  the  size  and  frequency  of  occurrence  of  earth¬ 
quakes  changing  markedly  early  on  15  June.  The  periodicity  of  the  large  events  during  the  13th 
and  14th  has  led  us  to  consider  elsewhere  various  deterministic  models,  including  tidal  trigger¬ 
ing,  in  discussion  of  the  cause  for  these  earthquakes.  Whatever  the  correct  physical  source 
model,  it  is  clear  that  Knopoff’ s  statistical  model  applies  in  this  case  where  we  found  that  the 
classical  frequency-magnitude  formula,  logN  =  a  +  bM,  did  not. 

J.  Filson 


E.  ENERGY-MOMENT  RATIOS  FROM  LONG-PERIOD  P  WAVES 

The  analysis  of  the  long-period  P  waves  from  a  Siberian  earthquake,  discussed  in  the  last 
1 

SATS,  has  been  extended  in  order  to  consider  the  complete  waveform.  First  motions  and  wave¬ 
form  elongation  had  previously  indicated  that  the  source  of  this  earthquake  was  a  vertical,  strike- 

slip  fault,  where  the  source  of  the  rupture  propagated  to  the  northwest  at  an  azimuth  of  325°. 

1 0 

Following  Haskell,  the  far  field  P-wave  displacement  (u(  )  from  a  linear  fault  with  unilateral 
rupture  is  represented  as 

U<*  =  Ra(r’0'  ^  MoT_1  tG(t)  "  G(t  _  T)]  (III-4) 

In  Eq.  (III-4) ,  R^  is  a  geometric  term  depending  on  the  epicentral  distance  r  and  the  orientation 
(0,  (p )  of  the  fault  plane  with  respect  to  the  recording  station.  Mq  represents  the  moment,  T  a 
characteristic  time  and  G(t)  the  time  variation  of  the  slip  across  the  fault.  In  Eq.  (III-4)  the 
characteristic  time  T  is  represented  by 

T  =  L  [C_1  -  o'1  cos  (©)]  (III—  5) 
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where  L  is  the  fault  length,  C  the  rupture  velocity  and  a  the  apparent  P-wave  phase  velocity. 

In  the  case  of  vertical  strike-slip  faulting,  0  is  the  azimuth  to  a  recording  station  measured  in 
the  horizontal  plane  from  the  direction  of  fault  propagation. 

Figure  III- 1 3  shows  how  one  might  gain  estimates  of  source  parameters  T,  Mq  and  energy 

E  in  an  idealized  case.  If  G(t)  is  a  step  function,  then  the  apparent  source  duration  would  be 

measurable  at  a  given  station  as  the  length  of  ground  displacement  u^  which  varies  as  G(t)  — 

G(t  —  T).  The  integral  of  u  should  be  proportional  to  the  moment  M  at  a  given  azimuth  and 
r  2  a  o 

the  integral  J  li  to  the  energy  in  the  waveform.  Ideally  these  integrals  should  vary  with  time, 

as  shown  in  Fig.  Ill- 1 3 .  The  product  of  the  ratio  E/Mq,  the  source  rigidity,  and  a  constant 

K(r,  0,  (p)  should  yield  the  apparent  effective  stress  acting  during  the  faulting  process.  Of  course 

we  cannot  expect  an  actual  seismic  source  to  give  rise  to  the  simplified  forms  of  these  functions 

shown  in  Fig.  Ill— 1 3 . 

Figure  III-14  shows  the  displacement  impulse  response  I (t)  of  the  long  period  WWSSN  instru¬ 
ment  at  BKS.  S(t)  is  the  trace  of  the  long-period  vertical  seismometer  at  BKS  [Berkeley]  and 
u(t)  is  the  ground  motion  gained  through  the  deconvolution  of  S(t)  with  I (t ) .  The  integral  of  u(t) 
with  time  is  shown  in  the  lower  left  of  Fig.  III-14.  The  long  term  level  of  this  integral,  if  cor¬ 
rected  for  radiation  and  transmission  effects,  is  proportional  to  Mq.  It  is  pointed  out  that  there 
is  a  base  line  problem  in  computing  the  integral  of  ground  displacement.  Out  of  some  30  wave¬ 
forms  to  which  this  technique  was  applied,  only  about  half  showed  a  clear  tendency  to  level  off 
at  some  constant  long  term  value. 

In  Fig.  Ill- 1 5,  H(t)  represents  the  response  of  the  instrument  to  an  impulse  in  ground  veloc¬ 
ity,  S(t)  is  the  same  seismogram  as  in  Fig.  III-14,  and  u(t),  the  ground  velocity,  is  gained  through 
the  deconvolution  of  S(t)  with  H(t).  The  time  variation  of  the  integral  /  u(t)2  at  stations  BKS  and 
NOR  is  shown  in  the  lower  left  of  Fig.  Ill— 1 5.  These  traces  are  quite  well  behaved,  usually  gain¬ 
ing  some  constant  level  within  20  sec  of  the  initial  motion.  The  NOR  and  BKS  integrals  are 
aligned  in  time  with  respect  to  the  first  detectable  P  motion.  The  more  rapid  rise  of  the  NOR 
integral  is  taken  to  be  the  effect  of  source  propagation  previously  discussed.  The  long  term 
level  of  these  integrals  is  proportional  to  the  energy  in  these  waveforms. 

Although  the  displacement  and  energy  integrals  discussed  above  are  yet  to  be  scaled  abso¬ 
lutely,  their  ratio  /  u(t)V /  u(t),  should  be  proportional  to  sin  (20)  in  the  case  of  the  Siberian 
earthquake  and  a  in  general.  In  order  to  test  the  former  hypothesis,  the  long  term  levels  of  the 
velocity /displacement  (or  energy/moment)  integrals  measured  at  20  sites  are  plotted  as  a  func¬ 
tion  of  0  in  Fig.  III-16.  The  solid  curves  in  this  figure  are  those  of  sin  (20),  drawn  to  approx¬ 
imate  the  azimuthal  variation  of  the  data  points.  To  be  consistent  with  the  fault  plane  solution, 
these  points  should  show  a  minimum  along  the  strike  of  the  fault,  325°.  This  is  not  the  case. 

The  points  are  nearly  a  maximum  at  that  azimuth  and  are  more  consistent  with  a  fault  plane  strik¬ 
ing  due  north  (or  due  west). 

In  an  attempt  to  resolve  this  contradiction,  the  amplitude  of  the  Love  wave  was  measured 
at  33  WWSSN  stations.  This  phase  was  well  dispersed  and  exhibited  a  clear  60-sec  component 
on  the  east-west  horizontal  seismometer  at  many  sites.  These  amplitudes,  corrected  for  spread¬ 
ing  and  attenuation,  are  plotted  in  Fig.  Ill — 1 7  and  should  vary  as  cos  (20)  with  azimuth.  As  the 
solid  curves  suggest,  these  amplitudes  have  a  maximum  near  due  north  and  west,  consistent  with 
the  integral  ratios  but  not  the  first  motion  data.  At  this  stage  in  the  study,  it  seems  that  the 
P-wave  integrals  and  the  Love-wave  amplitudes  are  giving  a  consistent  picture  of  the  total  fault¬ 
ing  process,  while  the  P-wave  first  motions,  in  this  case,  indicate  only  the  orientation  of  the 

initial  fault  breakage.  T 

&  J.  Filson 
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Fig.  Ill— 1 .  Fault  plane  solutions  of  Asia.  F  and  N  refer  to  solutions  of  Fitch^  and  Nowroozi.4  Base  map  is  a  seismicity 
map  published  by  the  National  Ocean  Survey  of  the  Department  of  Commerce.  Faults  are’main  deep  faults'  of  Yanshin.^ 


Fig.  III-2.  Slip  vectors  inferred  from  Fig.  Ill— 1 .  Base  map  same  as  Fig.  Ill— 1 . 
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Fig.  III-3.  Focal  sphere  for  Iranian  thrust  event.  Dilatations  are  open  circles, 
compressions  closed  circles. 
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Fig.  Ill— 4.  Relocated  earthquakes  for  Iranian  sequence.  Surface  breaks  are  also 
marked.  Open  circles  represent  less  reliable  locations.  M  is  the  main  shock. 
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Fig.  Ill- 5.  Short-period  P-waveforms  at  LASA  and  NORSAR  for  the  deep 
Hindu  Kush  event  of  20  January  1972.  The  expected  arrival  times  for  pP 
and  sP  are  marked.  The  records  should  be  read  like  Western  printing, 
left  to  right,  top  to  bottom. 
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Fig.  Ill— 6 .  Locations  of  the  20  January  1972  event  by  LASA,  NORSAR  and  NOS. 
The  array  locations  are  not,  of  course,  claimed  to  be  as  accurate  as  the  NOS 
location. 
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Fig.  III-7.  Short-period  P  at  College,  Alaska  (A  =  7  5°) 
for  the  events  of  2  9  January  1968  (upper)  and  20  January 
1972  (lower). 
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Fig.  III-8.  Long-period  records 
of  the  20  January  1972  event  at 
Hong  Kong  (HKC),  Istanbul  (1ST), 
Kevo  (KEV),  Nord  (NOR),  Baguio 
(BAG),  Chiengmai  (CHG)  and  Col¬ 
lege  (COL).  The  expected  arrival 
times  (approximate)  of  pP  and  sP 
are  marked  by  arrows. 
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Fig.  Ill- 9.  LASA  long-period  beam  to  the  20  January  1972  event. 
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Fig.  III-10.  Summary  of  the  seismicity  of  the  Galapagos  swarm.  Histogram 
at  top  shows  events  per  day  (Mg  >  3.0).  Middle  curve  is  seismic  energy  re¬ 
lease  on  a  logarithmic  scale.  Lower  line  is  seismic  energy  release  on  a  linear 
scale.  Note  linearity  of  energy  release  between  12  and  20  June. 


73 


Fig.  III-ll.  Results  of  numerical  evaluation  of  left-hand  side  (open  circles)  and 
right-hand  side  (solid  circles)  of  Eq.  (III-l)  over  the  period  13  to  19  June. 


DATE  JUNE  1968 


Fig.  III-12.  Details  of  seismicity  of  swarm  for  period  12  to  17  June.  It  is 
obvious  from  this  diagram  that  statistics  of  swarm  are  not  stationary. 
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118-2-10693(1)] 


a.  G(t)-G(t-T) 


d.  a  *  E  M0^ 

Fig.  Ill- 13.  Displacement  and  velocity  integrals  from  an  idealized  source, 
(a)  G(t)—  G(t  -  T)  represents  idealized  teleseismic  ground  motion  (uQ  )  from 
a  propagating  Heaviside  slip;  (b)  long-term  level  of  integral  of  this  ground 
motion  should  be  proportional  to  the  seismic  moment;  (c)  long-term  level 
of  integral  of  ground  velocity  squared  should  be  proportional  to  energy  in 
the  wave  form;  (d)  ratio  of  these  integrals  is  proportional  to  apparent  ef¬ 
fective  stress. 


Fig.  Ill- 14.  Computed  displacement  integral  at  BKS.  I(t)  is  displacement 
impulse  response  of  long-period  vertical  seismometer;  S(t)  is  P-wave  from 
a  Siberian  earthquake  written  by  that  instrument;  u(t)is  vertical  ground  mo¬ 
tion  gained  through  deconvolution  of  I(t)*u(t)  =  S(t). 


75 


Fig.  III-15.  Computed  velocity  integral  at  BKS.  H(t)  is  velocity  impulse 
response  of  long-period  vertical  seismometer;  S(t)  is  same  as  Fig.  Ill— 1 4 ; 
u (t )  is  vertical  ground  velocity  gained  through  deconvolution  H(t)*u(t)  =  S(t). 


Fig.  Ill- 1 6 .  Velocity-displacement  integral  ratios  as  a  function  of  azimuth. 
Open  dots  are  measured  ratios  of  velocity  to  displacement  integrals;  solid 
line  shows  theoretical  variation  of  this  ratio  with  azimuth. 
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Fig.  Ill— 1 7 .  Corrected  Love-wave  amplitudes  as  a  function  of  azimuth, 
dots  are  Love-wave  amplitudes  corrected  for  spreading  and  attenuation, 
line  shows  theoretical  variation  of  these  amplitudes  with  azimuth. 


Open 
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IV.  SIGNALS  AND  NOISE 


A.  SUPPRESSION  OF  LONG-PERIOD  NONPROPAGATING  NOISE 


The  work  on  the  prediction  and  suppression  of  long-period  nonpropagating  noise  described 

12  3 

in  the  last  two  SATS  *  has  now  been  completed  and  is  more  fully  described  elsewhere.  Out¬ 
lined  here  are  the  main  conclusions  of  this  study. 

Approximately  half  the  noise  observed  by  long-period  seismometers  at  LASA  is  nonpropagat- 
4 

ing;  that  is,  it  is  incoherent  over  distances  greater  than  a  few  kilometers.  However,  because 
it  is  often  strongly  coherent  with  microbarograph  data  recorded  at  the  same  site  at  the  same 
time,  a  large  proportion  of  it  can  be  predicted  by  convolving  the  microbarogram  with  some  trans¬ 
fer  function.  The  reduction  in  noise  level  using  this  technique  can  be  as  high  as  5  dB  on  the 
vertical  seismometer  and  higher  still  on  the  horizontals.  If  the  source  of  this  noise  observed  on 
the  vertical  seismometer  were  predominantly  buoyancy  in  the  instrument,  the  same  transfer 

function  would  be  expected  to  predict  the  noise  observed  in  the  output  of  the  same  microbarograph 

\ 

at  another  time  when  the  coherence  between  them  is  again  high.  Since  it  does  not,  one  is  forced 
to  conclude  that  the  noise  is  caused  by  atmospheric  deformation  of  the  ground. 

Since  so  much  of  this  noise  can  be  removed  merely  by  subtracting  that  component  which  is 
correlated  with  the  microbarogram  recorded  at  the  same  site,  the  deformations  must  be  caused 
largely  by  atmospheric  pressure  disturbances  in  the  near  vicinity  of  the  seismometer.  Any  far- 
field  atmospheric  loading  effects  must  be  quite  insignificant.  In  describing  the  loading  process, 
therefore,  we  would  expect  that  an  elastic  half-space  might  be  an  adequate  model  for  the  earth 
and,  since  the  noise  is  nonpropagating,  negligible  error  would  be  incurred  by  considering  the  de¬ 
formation  to  be  static.  Inspection  of  the  equations  for  the  displacement  due  to  a  two-dimensional 
sinusoidal  pressure  distribution  on  the  surface  of  a  half- space  shows  that  the  displacement  is 
directly  proportional  to  the  product  of  the  amplitude  of  the  pressure  distribution  and  its  wave¬ 
length.  If  the  pressure  distribution  were  stationary  and  moving  relative  to  the  ground  with  the 
wind  speed,  the  wavelength  of  the  noise  observed  by  the  seismometer  at  a  given  frequency  would 
be  directly  proportional  to  the  wind  speed.  One  can  see  why  the  transfer  function  which  works 

on  one  day  does  not  work  on  another:  the  wind  speeds  are  different, 

2 

A  hope  was  expressed  in  the  last  SATS  that  the  deconvolved  microbarograph  data  could  be 
weighted  with  the  deconvolved  anemometer  data  at  the  same  site  to  provide  a  time  series  from 
which  the  noise  could  be  predicted  with  a  time -invariant  transfer  function.  After  much  effort, 

3 

it  was  found  that  this  could  not  be  done.  The  wind  speed  at  any  time  is  generally  not  the  same 
over  space  even  within  the  coherence  distance  of  pressure  variations.  The  atmosphere  is  too 
turbulent.  In  other  words,  the  pressure  distribution  is  not  stationary  in  some  frame  of  reference 
which  moves  relative  to  the  ground  at  the  mean  wind  speed. 

Despite  the  turbulence,  or  rather,  because  of  it,  any  point  on  the  surface  of  a  half-space  is 
as  good  as  any  other  for  obtaining  time -averages  of  the  wind  speed  and  of  the  fluctuations  of 
microbarometric  pressure.  Using  the  half-space  model  and  applying  two  arbitrarily-chosen 
filters,  (1  —  e  yc)  and  e  r?c,  to  remove  any  constant  background  and  infinitely  high  wavenumbers, 

3 

where  c  is  the  wavenumber,  an  expression  has  been  derived  for  estimating  the  noise  power 
which  would  be  observed  by  the  vertical  seismometer: 
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P  (t)  is  the  mean-square  barometric  pressure 
2 

V  (t)  is  the  mean-square  wind  speed 

T  is  the  period  at  which  the  nonpropagating  noise  peaks  on 
the  vertical  seismometer 

v  is  Poisson's  ratio  for  the  half-space 

z  is  the  depth  into  the  half- space 

p  is  the  rigidity  of  the  half-space 

y  is  the  longest  observable  wavelength 

77  is  the  shortest  observable  wavelength. 


It  is  clear  that  the  logarithmic  term,  I  ,  becomes  infinite  if  y  is  infinite  or  if  77  and  z  are  both 

1 

zero.  In  other  words,  in  order  to  use  the  half-space  model  to  explain  band-limited  observations  of 
the  driving  forces,  the  pressure  disturbances  must  also  be  band-limited  in  wave-number  space. 

Equation  (IV— 1 )  was  tested  using  a  sample  of  data  lasting  about  hrs  beginning  at 
17h  20m  00s  on  23  August  1967  and  recorded  at  AO  at  LASA.  For  these  data,  77  and  y  were  esti¬ 
mated  in  the  following  way.  The  smallest  and  largest  periods  to  which  the  seismometer  can  re¬ 
spond  were  chosen  to  be  those  at  which  the  sensitivity  of  the  instrument  is  down  a  factor  of  10 
from  the  peak.  These  periods  are  about  9  and  80  sec  (see  Fig.  IV-1).  77  and  y  are  then  assumed 

to  be  equal  to  these  periods  multiplied  by  the  lowest  and  highest  wind  speeds,  5  and  10  msec  , 
recorded  during  this  time.  Thus,  77  =  45  m  and  y  =  800  m.  This  is  consistent  with  measure¬ 
ments  of  coherence  of  atmospheric  pressure  disturbances  made  by  Herron,  et  al?  From  Landers'^ 

measurements  of  P-  and  S-wave  velocities  in  the  sediments  of  UVSA,  Poisson's  ratio  v  is  esti- 

10  -2 

mated  to  be  about  0.3  and  the  rigidity  p  to  be  about  3.5  X  10  dyne  cm  .  From  the  microbaro- 

2  2  -4  2  2  -2 

graph  and  anemometer  outputs,  it  was  found  that  P  (t)  =  62.6  dyne  cm  and  V  (t)  =  37  m  sec 
From  the  coherence  spectrum  of  the  vertical  seismometer  and  the  microbarogram,  it  is  found 

that  T  =  40  sec.  Inserting  all  these  numbers  into  the  right-hand  side  of  Eq.  (IV-1)  and  setting 

2  18  2  2 
z  =  0,  we  find  E{u  }  =  723  x  10_1  m  .  From  the  vertical  seismogram,  we  find  that  u  = 

764  XdO’18™2. 

The  agreement  between  the  estimated  and  observed  noise  is  obviously  too  good  to  be  true. 

Uncertainties  exist  in  the  estimates  of  v  and  p^  which  may  affect  the  result  by  a  factor  of  2. 

The  factors  1^,  I^,  and  1^  are  somewhat  arbitrary  in  that  they  are  derived  from  arbitrary  as- 

4 

sumed  filters.  It  should  also  be  remembered  that  some  of  the  noise  is  propagating:  Capon 
showed  for  this  sample  of  data  that  the  nonpropagating  component  contributed  nearly  10  dB  more 
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to  the  noise  in  the  20-  to  4  0-sec  band  than  did  the  propagating  noise.  Therefore  only  about 
- 1 8  2 

70  x  10  m  of  the  observed  noise  can  be  propagating. 

It  is  interesting  from  a  practical  viewpoint  to  find  how  much  improvement  could  be  obtained 
by  burying  the  seismometer.  Figure  IV-2  shows  the  estimated  noise  power  E{u“}  plotted  as  a 
function  of  depth  z  using  the  data  discussed  above.  It  would  seem  from  this  example  that  the 
noise  power  of  this  particular  nonpropagating  noise,  generated  by  atmospheric  pressure  dis¬ 
turbances  with  wavelengths  on  the  order  of  45  to  800  m,  could  have  been  reduced  by  an  order  of 
magnitude  by  burying  the  seismometer  150  m. 

Sorrells,  et  al.,^  describe  a  very  careful  experiment  to  measure  noise  caused  by  local  at¬ 
mospheric  pressure  changes.  The  two  3-component  seismometers,  one  at  a  surface  site 
and  one  buried  183  m  in  a  salt  mine,  were  completely  isolated  from  any  measurable  buoyancy 
effects.  During  calm  periods,  the  noise  spectra  of  the  two  vertical  instruments  are  virtually 
identical.  During  a  windy  period  when  the  rms  pressure  fluctuation  was  12  jjbars  and  the  mean 

_  A 

wind  speed  was  7.8  msec  ,  the  noise  level  in  the  surface  instrument  increased  by  10  dB  in  the 
20-  to  100-sec  band,  while  that  in  the  mine  instrument  remained  the  same.  This  noise  condition 
is  approximately  the  same  as  discussed  above,  in  particular  the  observable  wavelengths  of  the 
pressure  distribution  are  about  the  same  —  perhaps  a  little  longer.  Figure  IV-2  predicts  that  at 
183  m  one  would  expect  to  get  a  little  over  11  dB  decrease  in  noise  level.  The  results  of  Sorrells, 

7 

et  al.,  thus  add  some  independent  confirmation  of  the  validity  of  Eq.  (IV- 1). 

There  would  have  been  some  future  in  predicting  and  eliminating  this  noise  had  it  been  pos¬ 
sible  to  obtain  a  time -invariant  transfer  function.  Since  this  appears  to  be  out  of  the  question, 
due  to  the  turbulence  of  the  fickle  wind,  there  remain  two  ways  of  suppressing  this  kind  of  noise. 
One  method  is  to  bury  the  seismometer  at  a  depth  of  about  150  m,  where  the  reduction  in  noise 
level  would  be  about  10  dB.  The  other  method  is  to  employ  an  array  of  seismometers  placed  no 
closer  together  than  1  km.  To  achieve  the  same  reduction  in  noise  level  as  would  be  obtained 
with  a  single  buried  instrument  would  require  10  surface  instruments.  At  LASA,  no  additional 
measures  need  be  taken  to  isolate  the  long-period  instruments  from  buoyancy  effects. 

A.  Ziolkowski 


B.  NOISE  AND  REGIONAL  RAYLEIGH  WAVES  AT  LASA 

The  question  addressed  was:  does  there  exist  a  minimum  in  the  long-period  noise  spectrum 
which  would  provide  an  optimum  band  for  detecting  Rayleigh  waves  from  sources  at  regional 
distances?  This  question  would  be  relevant  to  the  design  of  the  response  of  instruments  to  be 
used  in  a  detection  network  located  within  say  20°  of  a  seismic  region.  I  have  attempted  to  an¬ 
swer  this  question  by  comparing  the  average  long-period  vertical  noise  spectrum  recorded  at 
LASA  with  the  spectra  of  Rayleigh  waves  from  explosions  at  NTS  (A  =  12.5°). 

Noise  samples  of  140-sec  length  were  taken  at  sites  D1  and  F4  during  each  month  of  the  year 
1969.  The  amplitude  spectrum  of  each  noise  sample  was  computed  and  then  an  average  noise 
spectrum  taken  at  each  of  the  two  sites.  For  comparison,  the  amplitude  spectra  of  the  Rayleigh 
waves  from  the  NTS  explosion  PIPKIN  (m^  =  5.5)  were  also  computed  from  the  long-period  ver¬ 
tical  recordings  at  both  sites.  Both  the  signal  and  noise  spectra  were  computed  using  exactly 
the  same  processing  parameters,  such  as  data  length,  sampling  rate  and  taper  length.  The 
LASA  long-period  instrument  response  was  removed  from  all  of  the  spectra  in  order  to  give  a 
notion  of  relative  ground  motion  in  the  long-period  band. 
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The  results  of  this  analysis  are  shown  in  Figs.  IV-3  and  -4  (sites  F4  and  Dl,  respectively). 
Although  some  individual  noise  samples  were  bimodal  with  peaks  at  about  8  and  15  sec  this, 
apparently,  is  not  a  significant,  stationary  condition  at  these  sites  on  average.  The  average 
noise  sample  at  site  Dl  shows  a  slight  dip  at  0.08  Hz  (12  sec),  but  in  general  the  noise  field  at 
both  sites  increases  monotonically,  or  with  slightly  less  than  the  square  of  the  frequency  between 
0.04  and  0.15  Hz.  The  signal  spectra  at  both  sites  are  quite  similar  at  longer  periods,  increas¬ 
ing  as  the  first  power  of  frequency  in  the  band  0.02  to  0.1  Hz.  At  frequencies  higher  than  0.1  Hz, 
the  shapes  of  the  spectra  vary  significantly.  The  peak  signal  frequency  at  F4  is  about  0.9  Hz, 
while  at  Dl  it  appears  at  1.1  Hz. 

The  fact  that  the  peak  signal  frequency  at  one  site  occurs  at  a  dip  in  the  average  noise  field 
at  another  site  is  not  discounted,  indeed  on  certain  days  under  specific  noise  conditions  it  may 
be  quite  important.  On  average,  however,  the  advantages  of  shaping  the  instrument  recording 
response  of  a  regional  detection  network  expecting  an  inverse  relation  between  signal  and  noise 
near  a  period  of  12  sec,  do  not  appear  great. 

J.  Filson 
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Fig.  IV- 1.  Long-period  system  transfer  function. 


DEPTH  OF  BURIAL[z(m)] 


Fig.  IV- 2.  E{uz2}, 

the  following  data: 
T  =  40  sec;  jjl  =  3.5 


calculated  from  Eq.  (IV- 1),  as  a  function  of  z  using 
P2(t)  =  62.6  dyne2  cm  4;  V2(t)  =  37  m2sec~2; 

i0  -2 

x  10  dyne  cm  ;  v  -  0.3;  t)  =  45  m;  y  =  800  m. 
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Fig.  IV- 3.  Signal  and  average  noise  spectra  based  on  vertical, 
long-period  motion  at  LASA  site  F4. 


FREQUENCY (Hz) 


Fig.  IV-4.  Signal  and  average  noise  spectra  based  on  vertical, 
long-period  motion  at  LASA  site  Dl. 
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V.  DATA  HANDLING 


A.  INTERNATIONAL  SEISMIC  MONTH  PROGRESS 

Work  has  continued  on  the  International  Seismic  Month"'  during  the  past  six  months.  The 
1152  events  currently  on  the  bulletin  have  been  relocated  using  data  supplied  by  participants  and 
revised  epicenter  lists  distributed.  Participants  have  also  been  supplied  with  predicted  body- 
and  surface-wave  arrival  times  at  stations  of  particular  interest  to  them.  All  the  lists  and 
basic  data  for  these  1152  events,  including  all  arrival  times  used  to  verify  each  event,  have  been 
combined  and  entered  into  the  DADS  system  (see  Sec.  V-B).  Future  editing,  addition  of  new  data, 
calculation  of  new  derived  parameters,  and  listing  of  results  should  be  greatly  facilitated  by  this. 

Of  the  1152  events,  312  require  further  verification  to  be  certain  that  they  include  no  bogus 
events.  These  have  been  merged  with  an  additional  2000  potential  events  suggested  by  the  LASA 
and  NORSAR  detection  logs.  To  obtain  the  2000  detection  log  events,  some  automatic  sifting  of 
the  basic  logs  was  done  to  reduce  them  to  about  3800  entries  from  each  array.  Then  all  entries 
apparently  associated  with  one  of  our  840  verified  events  were  removed.  Also  removed  were 
obvious  multiple  detections  and  detections  with  insufficient  signal-to-noise  ratio.  A  serious  at¬ 
tempt  is  now  in  progress  to  verify  as  many  as  possible  of  the  more  than  2000  entries  on  this 
merged  list. 

A  large  body  of  data  has  been  assembled  and  more  will  be  obtained  to  try  to  enlarge  the 
number  of  verified  events.  All  reports  made  to  NOS  have  been  obtained  and  sifted  to  remove 
reports  associated  with  verified  events  and  to  remove  entries  from  stations  which  reported  fewer 
than  fifty  arrivals  during  the  month.  The  current  residual  in  this  list  is  about  12,000  arrivals. 
Similar  lists  of  several  thousand  as  yet  unassociated  arrivals  have  been  compiled  from  lists 
supplied  by  Canadian,  Swedish  and  British  colleagues.  Also  included  are  unassociated  arrivals 
at  selected  WWSSN  stations  obtained  by  picking  arrivals  using  film  records  rather  than  relying 
upon  lists  supplied  to  NOS. 

The  work  of  compiling,  integrating,  and  organizing  these  lists  of  unverified  detection  log 
entries  and  arrival  times  is  essentially  completed. 

Special  digitizing  and  processing  software  have  been  written  for  FM  analog  tapes  supplied 
to  us  by  the  United  Kingdom  Atomic  Energy  Authority.  The  objective  is  to  do  filtering  and  beam¬ 
forming  of  UK  array  data  to  aid  in  event  verification.  If  possible  and  profitable,  all  available 
Warramunga,  Australia  data  will  be  processed  with  these  programs.  Data  from  Gauribidanur, 
India  may  also  be  processed,  if  available.  Other  special  software  has  been  written  to  process 
event  tapes  from  LASA  and  NORSAR.  All  beams  and  subarray  beams  are  made  available  for 
visual  examination  with  and  without  filtering.  These  are  examined  selectively  as  questions  arise 
about  specific  events  or  detections.  ^  T  Lacoss 

R.  E.  Needham 
R.  M.  Sheppard 


B.  DATA  ANALYSIS  AND  DISPLAY  SYSTEM  (DADS) 

Development  of  a  new  interactive  software  system  to  display,  edit,  plot  and  otherwise  ma¬ 
nipulate  tabular  data  is  well  advanced.  Current  applications,  or  applications  under  development. 


*  Seismic  Discrimination  SATS,  Lincoln  Laboratory,  M.I.T.  (30  June  1972),  DDC  AD-748304. 
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include  management  of  ISM  data  and  of  our  entire  digital  tape  library  catalog.  It  is  also  being 
used  to  help  develop  NORSAR  station  corrections  for  our  computer  system. 

The  system  is  designed  to  deal  with  data  which  can  be  arranged  in  a  matrix  with  N  rows 
and  M  columns.  Size  limitations  are  NM^  400,000,  8,000,  and  500.  The  data  are 

stored  on  a  magnetic  tape  and  transferred  to  a  drum  at  the  start  of  each  session  with  the  data. 

All  changes,  additions  or  deletions  are  done  on  the  drum,  and  the  modified  data  base  is  stored 
on  tape  at  the  end  of  the  session. 

Each  column  of  the  matrix  is  assigned  a  name  by  the  user.  For  example  "DATE,"  "TIME," 

11 LAT,"  " LONG,"  and  "MB"  could  refer  to  columns  which  contain  date,  origin  time,  latitude, 
longitude  and  m^  values  for  seismic  events.  Each  row  of  the  matrix  could  represent  a  distinct 
event,  in  such  a  case,  and  entries  in  each  column  could  delineate  some  specific  seismic  param¬ 
eter.  Parameters  can  be  dates,  times,  integers,  floating  point  or  alphanumeric.  The  user  deals 
with  parameter  columns  by  their  names.  The  user  can  operate  with  all  columns  of  the  matrix 
or,  by  listing  names,  select  any  subset  of  columns  to  be  actively  manipulated  and  displayed 
during  a  session  with  the  data. 

Some  columns  of  the  matrix  are  denoted  selectors  and  serve  a  special  purpose.  Specifically, 
a  selector  defines  a  subset  of  the  rows  of  the  matrix  by  taking  on  a  special  value  (+1)  only  for 
rows  to  be  included  in  the  subset.  For  example,  a  selector  might  identify  all  events  in  a  list 
which  have  values  above  6.0.  Listing,  display  and  editing  can  all  be  done  subject  to  a  selec¬ 
tor,  so  that  extraneous  rows  do  not  obscure  the  process.  Selectors  can  currently  be  constructed 
row  by  row  manually,  or  by  specifying  rows  to  include  via  card  inputs.  Fortran  subroutines 
now  exist  to  allow  users  to  contruct  selectors  by  operating  on  other  parameters  and  a  general 
search  routine  will  be  added  shortly  to  expedite  this. 

Figure  V-l  shows  the  hierarchy  of  data  matrices  in  the  DADS  system.  A  large  matrix  of 
N  rows  and  M  columns  is  stored  on  our  drum.  An  active  submatrix  is  defined  by  a  list  of  param¬ 
eter  names  and  a  selector.  A  display  submatrix  of  the  active  submatrix  is  normally  visible  on 
a  CRT  display.  Finally,  one  element  of  the  display  submatrix  may  be  singled  out  and  put  into 
a  small  box.  The  box  can  be  moved  and  the  data  in  the  box  can  be  edited  at  will  from  a  teletype. 
The  location  of  the  display  submatrix  within  the  active  submatrix  can  be  moved  by  teletype  com¬ 
mands  or  via  analog  knobs  which  move  the  display  window  up,  down,  left  and  right. 

Figure  V-2  is  a  photograph  of  the  system  in  use.  In  the  figure,  the  basic  data  display  pro¬ 
gram,  which  is  the  heart  of  the  system,  is  being  used.  In  general,  all  other  programs  are 
called  from  the  display  mode  via  teletype  commands  and  they  return  to  the  display  mode  upon 
completion.  Currently  the  user  can  call  general  programs  to: 


(1)  Create  space  for  new  rows  or  columns  and  assign  names  to  the  columns. 

(2)  Enter  data  from  cards  into  any  active  submatrix  to  allow  batch  updates  or  additions 
to  the  data  base. 

(3)  Specify  a  set  of  columns  to  define  an  active  submatrix  by  entering  the  list  of  names 
from  cards. 


(4)  Define  a  subset  of  rows  via  a  selector  by  entering  a  list  from  cards. 

(5)  Print  any  active  submatrix  on  a  line  printer  in  a  standard  system  format. 

(6)  Output  any  active  submatrix  to  tape,  or  line  printer,  or  teletype  in  a  format  specified 
by  the  user. 
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(7)  Merge  two-ordered  list  into  a  single-ordered  list  by  switching  rows  in  the  data 
matrix. 

(8)  Generate  a  selector  which  defines  a  set  of  rows  as  the  union  of  other  subsets. 

(9)  Delete  rows  or  columns  from  the  data  matrix. 

In  addition  to  general  data  management  routines,  several  seismic  routines  have  been  incor¬ 
porated  into  the  system.  These  are: 

(1)  /ATVEL  to  generate  a  matrix  column  which  gives  event  arrival  times  at  a  fixed 
station. 

(2)  /ARRIV  to  obtain  body  wave  arrival  times. 

(3)  /DELAZ  to  calculate  distance  and  azimuth  to  a  station. 

(4)  /MB  to  calculate  body  wave  magnitudes  from  amplitude,  distance  and  period 

columns-  L.  J.  Turek 

A.  A.  Sholl 
R.  T.  Lacoss 
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Fig.  V-l.  Data  matrix  structure  of  Lincoln  DADS. 
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Fig.  V-2.  DADS  in  use. 
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VI.  GENERAL  SEISMOLOGY 


A.  OBSERVATIONS  OF  PcP  AND  P  PHASES  AT  LASA  AT  DISTANCES 

FROM  26°  TO  40° 

Steered  beams  of  short-period  P  and  PcP  phases  at  LASA  are  being  studied  to  determine 
the  nature  of  the  core- mantle  boundary  of  the  Earth.  In  spite  of  the  accumulation  of  vast  amounts 
of  short-period  data  at  single,  isolated  stations,  the  large  scatter  of  PcP  relative  to  P  ampli¬ 
tudes  prevents  a  clear  look  at  the  core- mantle  boundary. 

A  controversial  model  of  the  core- mantle  boundary  has  been  suggested  by  Buchbinder.  He 
reported  cases  of  reversed  polarity  of  PcP  phases  relative  to  P  at  single  stations  for  epicentral 
distances  less  than  32°,  where  PcP  passed  through  a  minimum.  Based  on  these  observations, 
Buchbinder  suggested  a  new  model  having  a  core-to-mantle  density  ratio  of  ~1.0. 

Here  we  report  our  observations  of  beamed  LASA  data  for  events  between  26°  and  40°,  and 

compare  them  to  theoretical  predictions  based  on  two  core-mantle  models,  one  given  by  Bullen 

2  1 
and  Haddon,  the  other  by  Buchbinder. 

The  results  of  the  analysis  of  forty-five  events  recorded  digitally  from  short-period,  vertical 
component  seismographs  at  LASA  are  reported.  Only  those  events  with  clean  P  and  PcP  arrivals, 
and  having  nearly  equal  periods  were  selected  for  study.  The  purpose  of  this  was  to  eliminate 
those  events  with  clearly  unequal  transmission  effects  or  radiation  patterns  along  the  P  and  PcP 
takeoff  directions. 

For  each  event,  steered  beams  across  the  entire  LASA  were  formed  for  both  P  and  PcP 
phases.  The  peak-to-trough  amplitudes  divided  by  period  (A/T)  were  calculated  for  each  beam. 
By  measuring  these  parameters  on  the  beams,  the  local  amplitude  variations^  caused  by  near 
surface  structure  below  LASA  are  minimized. 

Figure  VI-1  shows  the  P  and  PcP  phases  on  the  beams  for  4  typical  events.  The  similarity 
of  P  and  PcP  phases  is  quite  marked,  indicating  that  the  mantle- core  boundary  is  very  sharp 
where  sampled  by  our  data.  P  and  PcP  beams  for  a  selection  of  events  were  scaled  to  NOS 

Atyi 

magnitude  6.0  by  multiplying  each  amplitude  by  10  k,  where  Am^  equals  6.0  minus  the  NOS 
magnitude  of  the  event.  Figure  VI- 2  is  a  plot  of  the  normalized  P  amplitudes  (A/T)  in  m|a/sec 
as  a  function  of  epicentral  distances.  In  this  figure  and  all  others,  triangles  identify  events 
from  the  north,  and  circles,  those  from  the  south  of  LASA.  Assuming  a  spherical  radiation  pat¬ 
tern  at  the  source,  theoretical  P  and  PcF  amplitude  curves  were  calculated  for  earth  model 

4 

CIT208  of  Johnson.  The  amplitudes  were  smoothed  and  corrected  for  attenuation  at  1.0  sec  in 
the  mantle  using  the  short-period  Q  model  of  Archambeau,  Flinn  and  Lambert^  and  for  reflec¬ 
tion  at  the  free  surface.  The  P  amplitude  curve  normalized  to  give  A/T  values  for  NOS  magni¬ 
tude  6.0,  is  shown  as  a  solid  line  in  Fig.  VI-2.  The  observed  P  amplitudes  show  large  scatter, 
and  generally  are  larger  than  the  theoretical  P  amplitudes.  Lacoss^  observed  this  effect  for 
a  large  population  of  events  recorded  at  LASA. 

Figure  VI- 3  shows  the  observed,  as  well  as  theoretical  PcP  amplitudes  (A/T)  normalized 
by  the  same  scale  factors  used  for  the  P  amplitudes.  Curve  1  is  the  PcP  amplitude  (A/T)  cor¬ 
rected  for  geometric  spreading  and  mantle  attenuation  at  1.0- sec  period,  but  with  a  reflection 
coefficient  of  1.0  at  the  core-mantle  boundary.  The  effect  of  two  core-mantle  boundary  models 

is  shown  by  curves  B1  and  BK.  The  parameters  for  these  models  are  listed  in  Table  VI-1. 

2 

Curve  B1  is  the  PcP  amplitude  predicted  by  model  B1  of  Bullen  and  Haddon.  This  model,  with 
minor  modifications,  has  been  utilized  in  many  core  studies  as  a  standard  model.  Curve  BK  is 
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TABLE  VI-1 

CORE- MANTLE  MODELS 


Model 


BK-Buchbinder  (1968) 

mantle 

core 

Bl- Bullen  and  Haddon  (1967) 

mantle 

core 

BH1- Modified  Bl 

mantle 

core 

BH2-  Modified  Bl 

mantle 

core 

P  Velocity 
(km/sec) 

S  Velocity 
(km/sec) 

Core- to-  Manth 
Density  Ratio 

13.64 

7.29 

1 

7.51 

0.00 

13.64 

7.30 

1.74 

8.13 

0.00 

13.64 

7.80 

1.94 

8.63 

0.00 

13.64 

8.30 

2.14 

9.13 

0.00 

the  amplitude  predicted  by  the  core- mantle  boundary  of  Buchbinder.*  The  negative  sign  on 
curve  BK  for  A<  32°  indicates  a  negative  reflection  coefficient  at  the  core-mantle  boundary. 

The  observed  data  in  Fig.  VI- 3  are  widely  scattered  and  are  above  curve  Bl,  but  the  trend 
clearly  favors  model  Bl  over  model  BK.  The  slight  minimum  around  A  =  36.5°  cannot  be  con¬ 
sidered  reliable  in  view  of  the  observed  scatter  in  the  data.  We  neither  observe  any  minimum 
around  A  =  32°,  nor  the  change  in  polarity  below  32°,  as  predicted  by  Buchbinder. 

In  Fig.  VI- 4,  the  observed  and  theoretical  PcP/P  ratios  are  shown.  The  questionable  min¬ 
imum  at  ~36°  cannot  correspond  to  a  zero  reflection  coefficient  at  the  core-mantle  boundary, 
because  the  data  on  either  side  of  36°  have  the  same  polarity.  Except  for  the  five  indicated 
events,  the  data  show  no  change  of  polarity  or  amplitude  trend  with  distance  from  26°  to  40°.  A 

7 

similar  conclusion  was  reached  by  Bufe  in  a  study  of  PcP  phases  from  explosions  recorded  at 
WWSSN  stations. 

There  are  two  ways  to  obtain  theoretical  curves  for  PcP/P  ratios  which  are  slightly  higher 
than  curve  Bl.  One  is  to  modify  the  Bullen-Haddon  model  to  obtain  a  higher  reflection  coefficient 
at  the  core-mantle  boundary.  Two  examples,  BH1  and  BH2,  are  shown.  These  models  have  the 
parameters  given  in  Table  VI- 1.  Any  small  increase  in  mantle  shear  velocity,  core  density,  or 
core  compressional  velocity,  will  produce  an  increase  in  the  core-mantle  reflection  response. 
Curves  BH1  and  BH2  show  the  effect  of  increasing  all  three  of  these  parameters  simultaneously. 
Another  way  to  raise  the  PcP/P  ratio  is  to  increase  the  attenuation  of  P  relative  to  PcP.  The 

4 

Q  model  of  Archambeau,  et  al.,  used  for  the  theoretical  curves,  is  based  on  short-period  P 
data  from  explosions,  and  is  much  less  attenuative  than  models  based  on  long-period  data.  Be¬ 
tween  25°  and  40°  distance,  the  theoretical  curve  Bl  can  be  scaled  up  by  a  factor  of  2.0  by  using 

8  9 

more  attenuative  earth  models  such  as  those  suggested  by  Kanamori  and  Ibrahim  in  their  PcP 
studies. 

In  Fig.  VI- 5,  the  PcP/P  ratios  are  plotted  against  NOS  depth  of  focus.  The  wide  scatter  of 
data  above  100  km  probably  reflects  the  complexity  of  shallow  events,  as  well  as  uncertainties 
in  source  depth  and  attenuation.  Below  100  km,  the  data  scatter  appears  to  be  less,  but  this  is 
deceptive  because  the  density  of  points  is  also  less. 

In  conclusion,  our  observations  of  PcP/P  amplitude  ratios  at  LASA  between  A  =  26°  to  40° 
agree  with  theoretical  ratios  calculated  using  the  Bl  core- mantle  boundary  of  Bullen  and  Haddon 2 


94 


and  a  Q  structure  more  attenuative  than  the  short- period  Q  model  of  Archambeau,  et  al.  The 
wide  scatter  of  amplitudes  on  the  beamed  P  and  PcP  phases  prohibits  the  precise  determination 
of  the  core-mantle  boundary  parameters.  ^  ^  Frasier 

D.  K.  Chowdhury* 

B.  OBSERVATION  OF  FREE  OSCILLATIONS  OF  THE  EARTH 

1  0 

The  excitation  criterion  for  identification  of  spectral  peaks  of  normal  modes  of  the  earth 
enables  us  to  measure  the  eigenperiods  for  modes  which  could  not  be  identified  using  former 
techniques.  Figure  VI-6  shows  the  result  of  using  this  technique  to  identify  the  ^^q  mode,  ex¬ 
cited  by  the  Colombian  deep  shock  of  31  July  1970. 

An  exhaustive  analysis  of  the  data  available  for  this  shock  has  been  undertaken  and  approx¬ 
imately  100  new  modes  have  been  identified.  The  standard  deviation  of  the  eigenperiods  meas¬ 
ured  for  the  nS  modes  is  approximately  the  same  as  the  standard  error  of  the  mean  found  for 

the  Alaskan  earthquake  of  28  March  1964,  using  a  comparable  number  of  stations  (Eig.  VI- 7). 

1 2 

Derr's  eigenperiods  are  used  as  a  reference  in  Fig.  VI- 7. 

In  the  case  of  the  QSn  modes,  it  was  found  that  the  eigenperiods  are  systematically  smaller 
for  the  Colombian  earthquake  than  those  measured  for  the  Alaskan  event.  The  difference  is,  in 
general,  smaller  than  0.1  percent.  It  is  within  the  standard  deviation  of  the  measurements  at 
individual  stations,  but  larger  than  the  standard  error  of  the  mean  (s.e.m. ).  A  preliminary  con¬ 
clusion  of  this  result  is  that  the  lateral  heterogeneities  of  the  earth  are  such  that  for  a  given 
epicenter  and  the  actual  distribution  of  observing  stations  (the  WWSSN)  the  average  eigenperiod 
over  all  the  stations  depends  on  the  epicenter  location,  and  it  may  be  different  from  the  degener¬ 
ate  eigenfrequency  of  an  ideal  laterally  homogeneous  earth  by  an  amount  larger  than  the  s.e.m.  For 
the  JT  modes,  the  difference  is  about  3  times  larger  than  for  ~S  modes,  but  in  this  case  the 
eigenperiods  for  the  Colombian  shock  are  systematically  larger  than  those  found  for  the  Alaskan 
shock  (Fig.  VI-8). 

Those  results  clearly  imply  that  it  will  be  necessary  to  measure  eigenperiods  of  free  oscil¬ 
lations  excited  by  earthquakes  in  different  regions  of  the  world  in  order  to  estimate  the  degenerate 
eigenfrequency  of  the  ideal  laterally  homogeneous  earth,  and  in  order  to  study  the  spatial  distri- 
bution  of  the  lateral  heterogeneities.  j  Mendiguren 


*  Indiana  University  at  Fort  Wayne. 
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Fig.  VI-1.  Steered  LASA  beams  of  P  and  PcP  phases  for  4  events. 
Amplitudes  of  P  and  PcP  phases  have  been  equalized  to  demonstrate 
the  similarity  of  wave  shape  for  each  event. 
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DISTANCE  (deg) 

Fig.  VI-2.  Comparison  of  theoretical  and  observed  P  amplitudes  (A/ T), 
normalized  to  NOS  magnitude  6.0.  Theoretical  curve  (solid  line)  represents 
vertical  component  of  displacement /period  at  surface  after  corrections  for 
geometric  spreading,  attenuation  and  effect  of  free  surface. 
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Fig.  VI- 3.  Comparison  of  theoretical  and  observed  PcP  amplitudes  (A/T) 
normalized  to  NOS  magnitude  6.0.  Curve  1  is  for  unit  reflection  coefficient 
at  core-mantle  boundary,  curves  B1  and  BK  are  for  core-mantle  models 
of  Bullen  and  Haddon  (1967)  and  Buchbinder  (1968).  Negative  sign  on  BK 
curve  implies  a  reversed  polarity  of  PcP  from  P.  Solid  symbols  indicate 
observed  PcP  phases  with  opposite  polarity  to  P. 
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Fig.  VI- 4.  Comparison  of  theoretical  and  observed  PcP/P  amplitude  ratios. 
Curves  1,  Bl,  and  BK  are  for  the  same  models  described  in  Fig.  VI- 3.  Also 
shown  are  two  modified  Bl  models,  curves  BH1  and  BH2,  with  parameters 
given  in  the  text. 
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AMPLITUDE  RATIO  [pcP  (^)  /  P  (y)] 


Fig.  VI- 5.  Pc  P/P  ratios  plotted  against  depth  of  focus. 
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FREQUENCY  (cycles/sec) 

Fig.  VI- 6.  Result  of  application  of  excitation  criterion  to  identify  -S.q  mode 
on  colatitudinal  component.  ^  ^ 
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Fig.  VI- 7.  Fundamental  spheroidal  modes  qS_.  Difference  between  eigenperiods 
measured  for  Alaskan* *  and  Colombian  earthquakes  and  average  eigenperiods 
published  by  Derr.*  2  s.  D.  and  s.e.m.  indicate  standard  deviation  and  standard 
error  of  the  mean  found  for  Alaskan  earthquake  measurements. 
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RESIDUAL  (sec) 


Fig.  VI- 8.  Fundamental  torsional  modes  QTn  (see  Fig.  VI- 7  caption). 
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